de águas superficiais. Destaque aqui para a apli- 
cação da análise de variáveis regionalizadas (kriga- 
gem da piezometria e da transmissividade) e dos 
modelos matemáticos [9], [10, 1980]. 


6. ANÁLISE ECONÓMICA 


O planeamento do aproveitamento das águas sub- 
terrâneas, embora muito ligado nas regiões áridas a 
problemas de sobrevivência, difíceis de traduzir em ter- 
mos de economia, devem ser objecto de análise econó- 
mica. Os estudos de base nas diferentes escalas, as cap- 
tações, a exploração e conservação destas, e o controlo 
dos aquíferos, têm os seus custos. Em regime de eco- 
nomia de mercado é necessário medir esses custos e os 
benefícios resultantes da sua ocorrência, tendo em conta 
que o capital utilizado deve ser remunerado à taxa cor- 
rente das aplicações nas actividades a que as águas se 
destinam. 

Os princípios e métodos de análise aqui não têm 
nada de novo relativamente a outras actividades eco- 
nómicas, pelo que não serão descritos. Mencionam-se 
como exemplos de aplicação prática os casos de Cabo 
Verde [2, 1973], para ambiente rural (doméstico e agrí- 
cola), e da Península de Setúbal [10, 1980] para ambien- 
te urbano e industrial. 


7. LEGISLAÇÃO 


A necessidade de legislação como instrumento de 
coordenação das acções de planeamento e subsequente 
disciplina de utilização das águas subterrâneas é por 
demais evidente e justifica uma menção em capítulo 
especial, 

Não se pode, no âmbito deste trabalho, dar ao tema 
o desenvolvimento que ele merece, pois uma legislação 
adequada é tão importante para o êxito do planeamento 
quanto a sensibilização e empenhamento das popula- 
ções utentes da água subterrânea já mencionados em 
capitulo anterior. 

Essa legislação tem que ter em conta a relação entre 
as utilizações das águas superificiais e subterrâneas, en- 
tre os utentes do mesmo aquiífero e da mesma organiza- 
ção de conservação (em geral não será economicamente 
viável cada sistema de captação e elevação dispôr de 
organização própria de conservação) e entre os custos 
de controlo e protecção geral dos aquíferos e o preço 
da água ao consumidor. Deve permitir a tomada de 
medidas administrativas para a resolução dos proble- 
mas que as mencionadas relações implicam, e deve ain- 
da constituir instrumento adequado a dirimir questões 
entre as partes envolvidas no processo. 


8. CONCLUSÕES 


O Planeamento em águas subterrâneas é uma acti- 
vidade complexa e multifacetada. 

Do ponto de vista científico envolve conhecimento 
de diversas disciplinas como a meteorologia, a hidráu- 
lica, a química, a geologia, a matemática, a economia, 
etc. Do ponto de vista tecnológico, exige o domínio de 
capacidades que vão desde as mais rudimentares em 
captações até às mais sofisticadas, por exemplo, na uti- 
lização adequada de instrumentos matemáticos para mo- 
delagem de aquiferos. 

Contém uma forte componente socio-política e eco- 
nómica. Com efeito, ao destinar-se a fins que vão desde 
a sobrevivência das populações à actividade económica 
agro-pecuária e industrial, ao envolver extensões terri- 
toriais e investimentos muito significativos e ao exigir a 
coordenação com actividades afins, como as dos apro- 
veitamentos hidráulicos superficiais, não pode deixar 
de ser objecto de atenção, a todos os níveis, dos que 
têm responsabilidades administrativas e legislativas. E 
isto diz-se tanto em relação ao país onde a actividade 
de planeamento se desenvolve quanto em relação a orga- 
nismos internacionais [11, 1983]. 

Quem pois tem responsabilidades no planeamento 
de águas subterrâneas terá que atender à complexidade 
da sua natureza. Qualquer atitude miópica a este res- 
peito só pode concorrer para o fracasso de uma activi- 
dade que pode pôr ao serviço do homem recursos de 
valor muitas vezes insuspeitado pelos próprios que a ela 
se dedicam. 
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A CURVA LIMITE DE ESTAMPAGEM 
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RESUMO 


Com o presente trabalho pretende-se, através das 
CLE, analisar as condições em que actualmente se pro- 
duzem as bacias para lava-louças na Fábrica Portugal de 
modo a avaliar da possibilidade de se introduzirem alte- 
rações no fabrico que permitam a produção de bacias 
com paredes laterais ligeiramente cónicas por estampagem. 


ABSTRACT 


In the present work the production conditions of a 
sink in “Fábrica Portugal” are analised. 

Using the concept of the forming limit curves it is 
possible to access the modifications to introduce in the 
manufacturing process in order to produce a sink with 
slightly conical lateral walls by deep drawing. 


| - INTRODUÇÃO 


Os processos tecnológicos de deformação plástica de 
metais permitem obter peças a uma cadência elevada uti- 
lizando praticamente todo o material empregue. De entre 
esses processos, a ESTAMPAGEM, que permite a 
obtenção de uma peça de superfície não planificável a 
partir de uma chapa plana que manterá sensivelmente 
constante a sua espessura média, é um processo larga- 
mente utilizado. A estampagem é uma operação com- 
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plexa onde intervêm numerosos parâmetros, uns ligados 
à geometria da peça (cilíndrica, retangular, cónica, etc.), 
outros ligados às condições de operação (pressão do 
encostador, lubrificação, geometria e acabamento super- 
ficial da ferramenta, etc.) e outros ainda ligados às carac- 
terísticas do material seleccionado para a peça (coeficien- 
tes de encruamento e anisotropia, tensão de rotura, etc.). 

O encostador tem uma grande influência sobre o 
modo como se dá a deformação do material. Assim, 
quando ele exerce uma forte pressão sobre a chapa 
impedindo o seu escorregamento, o processo é designado 
por EXPANSÃO, sendo ambas as extensões maior e 
menor no plano da chapa positivas, enquanto que 
quando o encostador exerce uma menor pressão que, 
embora evite o engelhamento, permite o fluir do material 
sobre o punção, o processo é designado por RETRAC- 
ÇÃO, sendo neste caso a maior das extensões no plano 
da chapa positiva enquanto que a menor é negativa. 

Ao material estão associados problemas de instabili- 
dade durante o processo que vão provocar afastamentos 
em relação à peça pretendida motivados por defeitos de 
forma (engelhamento e recuperação elástica) ou pelo 
aparecimento de uma não uniformidade na distribuição 
da espessura ocasionada pelo aparecimento de um 
escoamento plástico localizado (estricção) que poderá 
conduzir a uma posterior rotura. Estes fenómenos limi- 
tam as possibilidades de deformação do material e são 
tanto mais relevantes quanto menor for a espessura da 
chapa e maior a sua tensão limite de elasticidade. A ten- 
dência actual é de utilizar cada vez mais materiais com 
estas características a fim de diminuir o peso dos compo- 
nentes e obter reduções significativas de consumos de 
energia e de matérias primas. Daqui se conclui da neces- 
sidade de uma análise profunda da influência daqueles 
fenómenos nos processos de deformação plástica. 

O conhecimento das deformações máximas que o 
material pode suportar pode ser obtido com o auxílio 
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das CURVAS LIMITE DE ESTAMPAGEM (CLE) (*). 
A CLE é uma representação gráfica no plano das exten- 
sões (e,, €;) (**) que o divide em duas zonas: uma zona 
inferior a que correspondem peças bem sucedidas e uma 
zona superior a que correspondem peças que “rompem” 
no decurso do processo de deformação. As CLE apare- 
ceram em 1965 com os trabalhos de Keeler (1) que as 
determinou em expansão biaxial, tendo sido completadas 
posteriormente por Goodwin (2) para o domínio da 
retracção. 

A CLE põe em evidência a não existência de um cri- 
tério simples de rotura ou estricção havendo pois uma 
função de e, e de e. 

O interesse prático das CLE assenta nas seguintes 
razões: 

— Nenhuma rotura ou estricção (***) se pode dar 
em pontos situados abaixo da curva. 


— Se por um processo de cálculo forem previstas as 
extensões e, e e, que terão lugar durante a estam- 
pagem então, de acordo com a posição desse 
ponto relativamente à CLE, determinar-se-á se a 
operação é possivel ou não. 

— Avaliação do grau de segurança com que se rea- 
liza a estampagem. 

— Qualificação dos diferentes materiais em função 
da sua aptidão para estampagem. 

— Análise da influência dos diversos parâmetros de 
estampagem. 


A CLE é traçada a partir dos resultados obtidos em 
diferentes ensaios laboratoriais de estampagem sobre 
chapas previamente marcadas com um determinado 
motivo: Rede. A precisão da CLE obtida depende do 
grau de precisão da marcação da rede e do tipo desta, do 
grau de precisão do processo de leitura, do tipo de 
ensaios realizados e, ainda, do próprio laboratório em 
que foram efectuados. 

Utilizam-se dois tipos de redes: as quadriculadas que 
são utilizadas para uma pequenissima dimensão do 
motivo (****), e mais frequentemente as redes formadas 
por círculos (afastados, tangentes ou entrelaçados) que 
dão de imediato a indicação das direcções principais no 
plano da chapa. Os processos mais usuais para a sua mar- 
cação na chapa são realizados com aresta cortante, com 
tinta ou por meios fotográficos, fotoquimicos, electro- 
químicos ou fotoelectroquimicos (6). 

De entre os diversos ensaios de estampagem (3), (4), 
(7). (8) e (9), os mais divulgados são o Bulge Test e 
ensaio Erichsen modificado, no domínio da expansão e o 
ensaio Swift no domínio da retracção. Existem outros 
ensaios que permitem determinar pontos da CLE em 
ambos os domínios, sendo de realçar de entre eles o 
ensaio Nakazima. 


A leitura das redes (4), (10) e (11) pode ser efectuada 
directamente ou por leitura num elemento intermédio 
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para o qual se “imprimiu” ou copiou a rede deformada. 
Para a obtenção de resultados significativos é neces- 
sário que: 


— A estricção ou rotura se tenha iniciado no circulo 
que se mede. 

— A estricção ou rotura passe pelo centro do circulo 
para que ele esteja disposto simetricamente, quer 
ao gradiente de deformação, quer à rotura. 


2. PARÂMETROS DAS CLE 


A caracterização do estado actual do conhecimento 
sobre as CLE é apresentado conjuntamente com os prin- 
cipais parâmetros que as influenciam. 


2.1. Dimensão da rede. Curva limite intrínseca 


Dos resultados experimentais (3), (4) e (5) as seguin- 
tes conclusões podem ser retidas: 


— A dimensão dos círculos depende da taxa de defor- 
mação local e não da dimensão da peça. 

— A CLE à estricção é pouco influenciada pelo diã- 
metro do circulo; a CLE à rotura depende bastante 
da dimensão do circulo. 

— Sanz e Grumbach (12) mostraram que o diâmetro 
do círculo influenciava muito o valor de e, e não 
influenciava o valor de e,. Foi assim possível 
obter as CLE intrinsecas (para um diâmetro de cir- 
culo nulo) que são independentes da trajectória de 
deformação (*****), 


(*) Na literatura Francesa: 
— COURBE LIMITE DE FORMAGE (CLF) 
COURBE LIMITE DEMBOUTISSAGE (CLE) 


Na literatura Inglesa: 


FORMABILITY LIMIT CURVE (FLC) 
— FORMING LIMIT DIAGRAM (FLD) 
FORMING LIMIT CURVE (FLC) 


(**) e, e e, são, respectivamente, as extensões principais, maior e 
menor, medidas no plano da chapa, a extensão segundo a espessura 
da chapa será designada por e,. As CLE são representadas quer 
em função das extensões nominais quer em função das extensões 
logaritmicas. 

(***) As CLE à estricção são mais dificeis de traçar devido à dificul- 
dade experimental em se determinar o início da estricção locali- 
zada, Para a sua determinação usam-se os métodos de Veerman, 
Bragard, Hecker e Kobayashi (3), (4), (5), conduzindo todos eles 
à mesma CLE à estricção, sendo a selecção do método a utilizar 
efectuada em função do equipamento disponível em cada 
laboratónio. 


(****) Os traços podem estar afastados de uma dimensão da ordem 
de 0,25 mm. 
(*****) Trajectória de deformação é a linha que representa a sucessão 
de deformações (e,, «,) que sofreu o circulo onde se produ 
ziu a estricção ou rotura final. 


2.2. Espessura da chapa 


A possibilidade de realizar com êxito uma operação 
através do aumento da espessura pode ser avaliada 
determinando o tipo de deformação predominante pois 
verifica-se que o aumento de espessura eleva a CLE à 
estricção significativamente no domínio de expansão (5), 
(13), (14) e (15). 


2.3. Trajectória de deformação 


As CLE referidas são obtidas em trajectórias simples, 
lineares, (3) e (5). Estas trajectórias lineares não represen- 
tam bem a situação real das estampagens industriais, 
particularmente quando se trata de estampagens em 
várias fases. As CLE determinadas em trajectórias com- 
plexas(*) dependem do modo como se deu a deformação 
na primeira fase (14), (15), (16) e (17). Da análise das 
CLE em trajectória complexa torna-se evidente a seguinte 
regra fundamental: começar sempre por uma deforma- 
ção em retracção antes de se fazer um passo em expan- 
são, para obter uma maior deformação global final. 

O CETIM(**) e o IRSID(***) apresentam uma 
explicação para esta influência através da análise das 
extensões equivalentes(****) (3), (4) e (5). 

A influência das trajectórias de deformação nas CLE 
convencionais é um dos principais problemas que se 
põem a uma maior utilização das CLE. 


2.4. Coeficientes de encruamento e de anisotropia 


O coeficiente de encruamento, n, traduz a capacidade 
do material em repartir as deformações sobre toda a 
superfície evitando, portanto, a localização da deforma- 
ção, melhorando desse modo a capacidade do material 
em suportar deformações por expansão, tanto mais 
quanto maior for o valor de n. Este facto justifica igual- 
mente a pouca influência de n sobre o coeficiente de 
estampagem, m(*****). 

Ao contrário do que acontece com o coeficiente de 
encruamento verificou-se uma acentuada influência do 
coeficiente de anisotropia, r(******), no valor do coefi- 


ciente de estampagem. Assim Sebaie e Mellor (18) verifi- 
caram que quanto maior for r, menor é o coeficiente de 
estampagem, o que é natural uma vez que o coeficiente 
de anisotropia traduz a capacidade de o material se opor 
ao adelgaçamento. 


Genericamente, pode-se dizer que quanto maior for o 
valor de r e de n, maior é a aptidão do material para 
estampagem, não alteram, sensivelmente, a posição da 
CLE mas sim as trajectórias de deformação (3), (4), (5), 
(19) e (20). Os resultados experimentais mostram que: 

— quando r aumenta a trajectória de deformação 
desloca-se para a esquerda (5) sendo esta influên- 
cia significativa em retracção e nula na zona de 
expansão. 


— quando n aumenta, a trajectória de deformação 
desloca-se para a direita e a rotura dá-se para 
valores de e, maiores, sendo esta influência signifi- 
cativa na zona de expansão e nula na zona de 
retracção. 

— As CLE em trajectórias complexas são tanto mais 
elevadas quanto maior for n (14) e (15). 

— A CLE em trajectória complexa sequência de 
tracção — expansão sobre quando r aumenta (14) 
e (15). 


3. TRABALHO EXPERIMENTAL 


A preparação das estampas para ensaio e a leitura e 
tratamento dos resultados foram efectuados no Labora- 
tório da Secção de Tecnologia Mecânica. As operações 
de estampagem decorreram na Fábrica Portugal utili- 
zando a prensa que produz os actuais lava-louças. 


3.1. Descrição do equipamento e da peça a obter 


A forma e as dimensões da bacia a estampar estão 
especificadas na Fig. 1. 

O material utilizado foi chapa de aço inoxidável 18/8 
(AISI 304) com as seguintes características: 

— composição química em % 


c< 0.08 Ps = 0.04 

81. =0 S max = 0.03 

Mn <20 Cr = 18.0/19.0 
Ni = 8.0/9.0 


— características mecânicas: 


or=53kg/mm” | Alongamento: 40% 

oe = 21 kg/mm” Dureza : 90 HRB 
— espessura: 

e = 0.7 mm 


Utilizou-se uma prensa hidráulica Miiller de duplo 
montante e duplo efeito com 500 toneladas no cabeçote 
principal e 350 toneladas no encostador. 


(*) Para simular uma trajectória complexa realizam-se dois 
ensaios de estampagem simples seguidos, ou seja, obtêm-se 
duas trajectórias lineares diferentes. 

(**) CETIM - Centre Technique des Industries Mécaniques, Sen- 
lis, France. 
(***) IRSID - Institut de Recherches de la Sidérurgie Française. 


s e 2 EST DE PCs TR=ro 
(****) Extensão equivalente € = da Vie 6 + lee? + tee? 


(*****) coeficiente de estampagem m: quociente entre o diâmetro do 


punção e o diâmetro máximo da estampa que pode ser estam- 
pada sem rotura por esse punção. Na literatura inglesa aparece 
frequentemente o Limit Drawing Rate (LDR) que é o inverso do 
coeficiente de estampagem. 


(******) Coeficiente de anisotropia r= E sendo: 


€ 
el - extensão logarítmica segundo a largura 
et - extensão logarítmica segundo a esnessura 
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Dos: Este desenho & aproximado pois, de 
facto, as paredes laterais têm um 
raio de curvatura de R=900 mm. 


FIG.l 


Bacia a estampar 


A ferramenta é formada por um punção e um encos- 
tador com pinos de centragem, montados na base da 
prensa e a matriz montada no prato superior da prensa, 
não havendo qualquer outro guiamento entre punção e 
matriz que não seja o resultante do próprio guiamento 
da prensa. O punção e o encostador são integralmente 
em aço enquanto que a matriz tem um anel em bronze 
alumínio para diminuir os atritos e os desgastes. 

As características geométricas determinantes no pro- 
cesso de estampagem são as seguintes: 


— raio do canto do cunho r,. = 20 mm 
— raio do canto da matriz r., = 2 mm 
— folga na zona dos cantos j = 0.99 mm 
— folga na zona dos lados j = 0.88 mm 


3.2. Ensaios 


Realizou-se um primeiro ensaio experimental para 
definir o processo de marcação das redes, a dimensão 
dos círculos, o processo de leitura e as secções típicas a 
estudar. Os resultados deste ensaio conduziram a adop- 
tar marcações de círculos de 10 e 15 mm de diâmetro 
com canetas de acetatos e segundo as secções típicas 
indicadas na Fig. 2. A leitura das elipses foi efectuada 
com uma luneta micrométrica equipada com uma escala 
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graduada e ampliação (7X) à excepção das que se encon- 
travam na zona dos cantos do “cunho” e “matriz” que 
foram copiadas para uma fita adesiva transparente e pos- 
teriormente lidas com a mesma luneta micrométrica. 

Embora a bacia tenha apenas um eixo de simetria, 
considerou-se a existência de quatro quadrantes (ITQ, 
HQ, HIQ e IVQ) havendo em cada um deles duas sec- 
ções típicas “simétricas” à excepção, obviamente, do can- 
to/secção central. 


== Contorno teórico da estampa plana 
Contorno real da estampa plana 


— ..—  Gacções estudadas 


pç INTERMÉDIA 
10, lado IIQ d=10 (15) 
[EMO TRANSIÇÃO 
10, lado ILO d=10 
(CANTO/SECÇÃO INTERMÉDIA 
So Ty lado TQ d=10 
e Mira ORNT SEER CENTRAL 
CANTO /SECÇÃO INTERMEDIA 
1Q lado IVO d=10 
ANTO-LADO/SECÇÃO TRANSIÇÃO 
IQ, lado IVO d=10 
f LADO/SECÇÃO INTERMÉDIA 
IQ, lado IVO d=1D (15) 


IO 


* Os círculos de 10 mm de 
diâmetro são, eventual- 
mente, entrelaçados. 


FIG.2 


Geometria da estampa plana. Rede e nomenclatura 


utilizadas nos ensaios nº 2, 4 e 5 


Após o ensaio experimental (EN. 1) realizaram-se 
mais quatro ensaios nas seguintes condições: 

Ensaio 2 (EN.2): utilizou-se uma rede formada por 
circulos entrelaçados no canto/secção central e canto/ 
secção intermédia(*), estando as condições do ensaio 
especificadas na Fig. 3. Utilizou-se uma pressão elevada 
no encostador(**) e sem lubrificação de modo a promo- 


(*) cerca de 600 circulos por ensaio. 
(**) As condições industriais de estampagem são: (as do EN.1) 
Força no encostador: 200 ton. 
Força máxima de estampagem: 60 a 80 ton. 
Lubrificação: tal como nos ensaios EN.4 e EN.S. 


ver o aparecimento de rotura e assim obter pontos de 
rotura no plano das extensões (e,, €,). 

Ensaio 3 (EN.3): ensaio idêntico ao anterior em que 
se diminui, ligeiramente, a pressão do encostador e 
lubrificou-se apenas a chapa do lado da matriz e somente 
nas zonas que irão formar os cantos da bacia. Tal como 
no anterior houve rotura. 

Ensaio 4 (EN.4): utilizou-se uma rede formada por 
círculos tangentes, sendo as condições do ensaio as espe- 
cificadas na Fig. 4. Neste ensaio simulou-se a pressão 
mínima do encostador sem engelhamento. 

Ensaio 5 (EN.5): utilizou-se uma rede formada por 
circulos entrelaçados no canto/secção central e canto/ 
secção intermédia, sendo as condições de ensaio as espe- 
cificadas na Fig. 5. Com este ensaio simulou-se o 
aumento de pressão do encostador sem provocar rotura. 


— rotura | Velocidade média 8.41 mm/s 


Eh pressão encos tador “*de estampagem: 0.5 m/min 


O posicioncdores Força encostador: 250 ton, 


Força máxima estampagem: 50 ton. 
Sem lubrificação. devido à rotura não houve a for- = | pressão ercostador 
ça máxima para vincar o fundo. “O positionadores 


FIG.3 
Condições do ensaio 2 (EN.2) 


= | Pressão encostaodor 
O posicionadores 


| Velocidade média 7.48 mm/s 

> » de estampagem: 0.67 m/min 

SS, al 
CREA] Força encostador: 230 ton, 
Força maxima estampagem: BO ton. 
1BO ton. no fim para vincar fundo 
À Lubrificação: Castrol iloform 
PS 158 

lado matriz: toda a chapa 
lado punção: só na região (ea 


"* Velocidade média 11.21 m/s 


de estampagem: 0.67? m/min, FIG.5 
Força encostador: 150 ton. Condições do ensaio 5 (EN.5) 
Força máxima estampagem:70 ton. 
250 ton. no fim para vincar Fundo 
Lubrificação: Castrol ilaform 
PS 158 : : o : 
é 4 Ledo petrias todo é capa Estes ensaios permitem analisar o efeito da pressão 
E L ção ] ã , . pe 
Sses KO] tado punçãos só na região (Ga do encostador (EN.4 e EN.5), da lubrificação (EN.2 e, 
EN.5) e ainda, do modo como a pressão se distribui pela 
FIG.4 chapa num mesmo ensaio. | 
ondições do ensaio EN. | ] adas var s signifi- 
Cond d 4 (EN.4) Neste estudo não foram detectad ações signifi 


89 


cativas na posição das frentes de deformação, F.D.(*), 
entre o EN.2 e o EN.3 pelo que não se apresentam os 
resultados obtidos no EN.3. 


4. RESULTADOS 


A influência da pressão do encostador manifestou-se 
sob dois aspectos: 

a) Uniformidade ou irregularidade na distribuição de 
uma dada pressão do encostador nas diversas 
regiões de uma estampa plana. 

b) Valor, maior ou menor, da pressão do encostador. 


F.D. Canto/sécção central 


IQ e IVO EN,4 


F1G.6 


As figuras 6 e 7 e as figuras 8 e 9 revelam o anda- 
mento idêntico entre as F.D. dos IQ e IVQ entre si e dos 
HQ e HIQ entre si. As F.D. do IQ e IVQ apresentam em 
relação às F.D. do IIQ e IIIQ uma maior expansão na 
zona do fundo da bacia e na parte cilíndrica, retracção, 
evidencia-se um aumento da extensão e, e uma ligeira 
diminuição (em módulo) da extensão e,. Estas evoluções 
resultam da assimetria com que a pressão exercida pelo 
encostador se distribui na chapa. Verificou-se a existên- 
cia de pressão mais elevada no IQ e no IVQ. 

O IVQ e o IHQ apresentam, na zona de retracção, 
um e, mais elevado do que o IQ e o IIQ, respectiva- 
mente, uma vez que a bacia é ligeiramente mais alta do 
lado do HIQ e IVQ. 


(*) Neste trabalho designar-se-á por FRENTE DE DEFORMAÇÃO 
(F.D.) a linha que melhor se ajusta aos pontos representativos, no 
plano das extensões (e,, €,), das deformaçães sofridas pelos diferen- 
tes pontos do material ao longo de uma dada secção da bacia. Estas 
extensões (e,, €,) foram medidas no fim da operação de estampagem. 
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E, 


——— 1TIQ 
== TI1Q 


F.D. Canto/secção central 


HQ e TIIO EN.4 


Mesmo sem analisar os resultados do EN.4 e EN.5 
entre si, pode-se avaliar da influência da pressão do 
encostador analisando, num único ensaio, secções idênti- 
cas mas sujeitas a pressões diferentes. 

A análise anterior vem reforçada pela análise compa- 
rativa das figuras 6 e 8 e figuras 7 e 9, onde se evidencia 
para o caso do EN.5 (maior pressão do encostador): 


— aumento da expansão no fundo 

— tendência ao aumento do e, na parte cilindrica, 
retracção 

— ligeira diminuição, em módulo, do e, na mesma 
parte cilíndrica. 


F.D. Canto/secção central 


D.? 
IQ e IVO EN.5 


F.D. Canto/secçao central 


- e 


IIQ e IIIQ EN.5 


FIG.9 


É evidente a melhor concordância entre as F.D. do 
EN.5 do que a que existe para as F.D. do EN.4 que se 
explica pela utilização de círculos entrelaçados no EN.5. 

A análise da influência da pressão do encostador 
atrás exposta verifica-se igualmente nas restantes secções 
(6), havendo uma menor concordância entre as F.D. do 
que a que existe no canto/secção central, uma vez que 
sendo menores as deformações vem ampliado o efeito da 
precisão de leitura e de marcação das redes, assim como 
da assimetria das secções em estudo. 

O factor de ampliação das deformações, Quadro 1, 
não apresenta variações significativas. 

A influência da não uniformidade da distribuição da 
pressão é também evidente pela análise dos resultados do 
EN.2 que se apresentam em seguida. 


QUADRO 1 
FACTOR AMPLIAÇÃO EN.S 


(*) 
SECÇÃO 


Eantoiidssgão central feto 5.8 


|Canto/Secção inter- 


né pes Jus | 1.6 Jus 
media 
Canto-lado/Secção | 
pair Vol Tal 1.8 
de transição | 
Lado/Secção inter- 
| média ' 4 (A 


(*) pErssdsrdo sa os valores de E| e E» 
do lado/secção intermédia como va- 
lores de referência. 


A influência da lubrificação é avaliada pela análise 
dos resultados do EN.2, figuras 10 e 11, comparativa- 
mente aos do EN.5, tendo-se observado: 


— Significativo aumento da expansão no fundo. 

— Grande aumento da extensão e, na região dos e, 
= (), resultando deste facto o aparecimento de 
rotura. 


FIG.10 


F.D. Canto/secçao central IQ e IVQ EN.2 


e, 
7 D.5 
EF ER 
: A 0.4 
de a 
ÇA “Sa 
a Ea 
" D. Md 
Ma ” 
+ o 
“ : 
“e 0.2 E: 
— IIQ e 
sm mim TITO | ; 
D.1- 
Ea 
=.3 =D. =.1 O O 1 De? 
FIG.11 


F.D. Canto/secção central IIQ e IIIQ EN.2 
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— Os pontos de rotura, obtidos no EN.2 para o can- 


to /secção central e canto / secção intermédia locali- 
zam-se na região dos e, = 0, Fig. 12. 


— Nos quadrantes onde a pressão era mais elevada, 


IQ e IVQ, apareceram roturas, que pela análise 
dos resultados se conclui surgirem primeiro no IQ 
e depois no IVQ tendo-se propagado e vindo a 
encontrar-se. 


— À rotura deu-se para uma penetração do punção 


superior a r,. + Fm que se pode explicar pelo 
facto de se estar a utilizar um aço inoxidável aus- 
tenítico que encrua rapidamente quando defor- 


tria na distribuição da pressão mostrou ter bas- 
tante importância, nomeadamente para determi- 
nar as secções onde se inicia a rotura, pelo que há 
que ter cuidado especial na manutenção das 
prensas. 


— O aumento da pressão do encostador provoca: 


a) aumento da expansão no fundo da bacia; 

b) aumento de e, na parte cilindrica (retracção); 

c) diminuição do módulo de e, na parte 
cilindrica; 

d) quando a pressão aumenta, pelo facto de o 
ejaumentar na região dos e, = 0 há a possi- 


bilidade do aparecimento de rotura; 
e) não se notaram alterações significativas no 
valor do factor de ampliação das deforma- 


mado plasticamente a frio. 


E, a 
ções. 
0.5 — A não utilização de lubrificação mostrou-se parti- 
2 cularmente inadequada, pois: 
ee a) aumenta muito a expansão do fundo da 
Ke : A R bacia; | | 
b) aumentam muito as extensões e, na região 
dos «e = 0 com o consequente perigo de 
0.5 rotura; 
"Pontos de rotura" obtidos c) aumenta o desgaste da ferramenta; 
0.2 d) aumenta a temperatura de trabalho. 


no ensaio EN.2 no canto/ — Atendendo à pequena espessura da chapa e aos 


ensaios realizados não parece possível, mantendo 

0.1 as características geométricas fundamentais do 
punção e da matriz (r.. = 20 mmer,, = 2 mm),a 
inclinação da parede lateral da peça pois o 
aumento da pressão do encostador necessário 
para evitar o engelhamento provocaria o apareci- 
mento de rotura. 

— Aconselha-se a mudança da geometria da bacia 
no sentido de aumentar o raio do canto da matriz 
(rm = 5 mm) e diminuir o raio do canto do 
cunho (r.. = 4 mm), bem como aumentar o 
raio de curvatura das partes laterais da bacia. 
Com esta alteração geométrica e um ligeiro 
aumento da espessura da chapa, os resultados 
experimentais levam a concluir ser possível obter 
peças com paredes laterais ligeiramente cónicas 
numa única estampagem. 


secção central e canto/ 


secção intermédia 


FIG. 12 


5. CONCLUSÕES 


A peça em estudo apresenta um único eixo de sime- 
tria o que criou dificuldades adicionais em virtude das 
alterações que produz nas F.D. a não simetria biaxial. 

Como principais conclusões a reter tem-se: 

— Um raio do canto da matriz, r.,., excessivamente 

pequeno provoca um grande aumento das exten- 

sões €,; 

— O r., muito pequeno associado a um grande raio 
do canto do cunho, r., provoca um grande 
aumento das tensões para que se inicie a estampa- 
gem da aba, com o consequente aumento da 
expansão na zona do fundo da bacia; 

— Necessidade de se utilizarem ferramentas bem 
guiadas. Da não existência de guiamento nesta 
ferramenta resultou uma eventual não uniformi- 
dade na distribuição da folga e na “centragem da 
chapa”, bem como uma assimetria na distribuição 
da pressão que o encostador exerce. Esta assime- 
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THERMAL NOISE IN SINGLE INJECTION 
SOLID STATE CYLINDRICAL DIODE WITH 
THERMAL FREE CARRIERS 


S. A. HUSAIN* 
Y.K. SHARMA* 
N.Kk. SHARMA* 


The analytical expressions for the thermal noise 
have been evaluated for the single injection current 
flow in solid state cylindrical diode with thermal free 
carriers. The regional approximation and salami 
methods are applied to obtain the results. It is shown 
that the thermal noise is suppressed by the space 
charge. 


|. INTRODUCTION 


The studies on the electrical conduction and noise 
behaviour for different current flow geometries have 
been made from a long time [1-5]. The presence of 
various physical parameters plays an important role 
on the space-charge-limited current flow in insula- 
tors. The scattering of current carriers with the 
lartice imperfections gives the thermal noise across 
the diode. The thermal noise is generated in the 
complete range of current-voltage characteristic which 
is modulated in the presence of various physical para- 
meters. The geometry of the diode is important for 
the magnitude of the current flow. Therefore, the 
thermal noise has been estimated in the cylindrical 
diode with thermal free carriers. The effects of space 
charge on the current injection and thermal noise 
have been studied with the help of regional appro- 
ximation and salami methods, respectively. It may be 
noted that the thermal free carriers contribute to 
the current flow but not to the space charge. 


* Department of Physics, R.B.S, College, Agra, India. 
Original recebido para publicação em 5/3/85 
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2. CURRENT INJECTION 


Let us consider a single injection solid state cylin- 
drical diode which contains a significant density of 
thermal free carriers n,. The general equations for 
the current flow and Poisson's law for the problem 
are given by 


I=2reurn(r) E(r) (1) 
E 1 d(rE) = 
Er E do MD-n (2) 


where I is the current per unit length of the cylinder, 
e the magnitude of the electronic charge, wu the free 
electron drift mobility, E the electric field intensity, 
n(r) the concentration of free carriers at radial dis 
tance r, e the permittivity of the insulator and r 
the radial distance inside the insulator varying from 
cathode radius r. to anode radius r.. 

The insulator is divided into two separate regions 
with the help of regional approximation me- 
thod [1-3]. The region I is space-charge-limited and 
the region II is Ohmic. These regions are well se- 
parated by the imaginary transition plane r,. The 
general equations for both the regions are obtained as 


Region I (r.<r< rn: 


I = 2x e urn(r) E(r) (3) 
en 1 drB) — 
Ea] de a 


Region II (r,<Ir< ff): 


I=2revuvnrE(r) (5) 
[E A das) À O (6) 
- Tr dr 


where r, is the transition plane at which the thermal 
free carriers are equal to the injected carriers. The 
voltage applied across the insulator is given by 


is f E(r) dr (7) 
Lc 


3. THERMAL NOISE IN THE DIODE 


The thermal noise is generated in the complete 
range of current-voltage characteristic. The expres- 
sions are evaluated with the help of results obtained 
for the current injection. The salami method of noise 
calculation [4,5] has been applied in this section. 
The whole insulator is divided into small slabs by 
the concentric cylinders. The noise source for ther- 
mal noise is present in each and every section of the 
diode. The thermal noise is estimated in one slab of 
the insulator. The integration from cathode radius 
to anode radius gives the total amount of thermal 
noise present in the insulator. 

According to Van der Ziel [4,5], the mean 
square noise voltage in a frequency interval Af is 
given by 


vi =4kTR AÍ (8) 


where 
R= = (9) 


In the above equations, k is the Boltzmann's cons- 
tant, T the lattice temperature and R the D.C. 
resistance of the diode. The thermal noise is estimat- 
ed in the complete range of current-voltage characte- 
ristic as given in order: 


a) — Low Injection Level: 


The injection level is very low, therefore the 
region II is only present in the insulator. The current- 
“voltage characteristic of this regime is derived from 
equations (5) and (7) as 


I=2reva, | In ( E )| EE (10) 


The noise resistance of the diode at low injection 
level is obtained from equations (9) and (10) 


R=——————— (11) 


From equations (8) and (11), the thermal noise in 
the cylindrical diode at low injection level is given by 


==. M (12) 


b) — Medium Injection Level: 


Both the regions I and II are present in the 
insulator. The current-voltage characteristic of the 
diode at medium injection level is derived from the 
equations (3) — (7) 


I 21u N %4 17 
+ so | In | eng, Lar] | (13) 


which directly gives the noise resistance of the diode 
as 


= 1 à Hr 3 | 1 [ 
"eum 2el 2'] 2reul 


a. E Et MD 
este [fo ( o 


The mean square noise voltage for this transition 
regime is given by equation (8) where R follows 
from equation (14). 


c) — High Injection Level: 


Only the space charge region I is present inisde 
the insulator. The current-voltage characteristic vf 
pure space charge current is derived from equa- 
tions (3), (4) and (7) as 


| o V 2 
I=2reu To. 4 E 3 (15) 
(a) = 
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The noise resistance of the diode is evaluated from 
equation (15) as 


ss EfO A e] tá 
R= | (=) gt | 


The mean square noise voltage for space charge re- 
gime may be evaluated from equations (8) and (16) 


2kT Af E (E 3 2 
=" lr+t—e([—— )Jj=-——+€ 
pe MET (E) 3] am 


which reveals that the thermal noise is inversely 
proportional to the applied voltage in the space- 
-charge-limited current flow in single injection solid 
state cylindrical diode with thermal free carriers. 


4. DISCUSSION AND CONCLUSIONS 


In the present analysis, the expressions for ther- 
mal noise have been evaluated for the three current- 
-voltage regimes of single injection solid state cylin- 
drical diode with thermal free carriers. The 
expression of thermal noise at low injection level 
shows that the noise is a constant quantity in Ohmic 
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regime where the current is directly proportional to 
voltage. The effect of thermal free carriers appears 
in the medium injection level of current where both 
the regions are present in the insulator. The thermal 
noise decreases from the higher value which is obtain- 
ed at low injection level. The thermal noise de- 
creases rapidly at high injection level because the 
thermal noise is inversely proportional to the applied 
voltage in space charge regime. The thermal noise 
is suppressed by the injected carriers because the 
concentration of thermal free carriers is compensat- 
ed by the injected current carriers. The important 
result of the present work is that the space charge 
suppresses the thermal noise. 
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TRANSITIONS BETWEEN EQUILIBRIUM 
CONFIGURATIONS IN A SOLID- LIQUID 


SYs TEM 


PAULO CARREIRA * 
M. AMARAL FORTES** 


RESUMO 


Ás possíveis configurações de equilíbrio de um 
sistema compreendendo corpos sólidos e um líquido são 
classificadas em termos do número e tipo de regiões 
líquidas. Quando se atingem determinadas configu- 
rações críticas, ocorrem transições no sistema, as 
quais originam diferentes regiões líquidas. As tran- 
sições podem resultar quer do contacto de regiões 
liguidas (entre si, ou com os sólidos) quer de «ultra- 
passagem» do limite de estabilidade de uma dada 
região. Com base em argumentos termodinâmicos e 
dinâmicos simples, prevê-se a configuração final no 
caso de um sistema com simetria cilíndrica, na ausên- 
cia da gravidade. Nas transições por contacto, a con- 
figuração final depende fundamentalmente da curva- 
tura das superfícies líquidas que contactam; e os 
estados intermédios na transição, considerados como 
estados de equilíbrio, têm valores decrescentes da 
energia de Helmholtz. 

O grau de conexão do líquido diminui, em regra, 
nas transições por contacto mas pode aumentar nas 
transições por «ultrapassagem» do limite de estabili- 
dade. 


ABSTRACT 


The equilibrium configurations of a system com- 
prising solid bodies and a liquid are classified in 
terms of the number and type of the liquid regions. 
At particular critical configurations, transitions occur 
that lead to different liquid regions. The transitions 
may be produced by contact of liquid regions or by 
overshoot of the stability limits of a region. The 
final configuration is predicted from simple dynamic 
and thermodynamic arguments for a system with 
cylindrical symmetry, in the absence of gravity. In 
contact transitions the final configuration depends 
essentially on the curvature of the fluid interfaces 
in contact and is reached through a path of decreasing 
equilibrium Helmholtz energy. The degree of connec- 
tion of the liquid generally decreases in contact trans- 
itions but may increase in overshoot transitions. 


| — INTRODUCTION 


A closed system comprising given volumes of two 
immiscible fluids and solids of a given geometry may 


* Laboratório Nacional de Engenharia e Tecnologia Industrial (LNETI) — Paço do Lumiar, Lisboa. 
** Departamento de Metalurgia, Instituto Superior Técnico; Centro de Mecânica e Materiais da Universidade Técnica de Lisboa (CEMUL). 


Av. Rovisco Pais, 1000 Lisboa. 
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L 


Fig. 1— Types of liquid regions in one-rod (a) and two- 
«rod (b) systems. Non-symmetric regions are examplified 
in (c) 


admit several alternative equilibrium configurations, 
each corresponding to a particular distribution of the 
two fluids. It will be assumed that one of the fluids 
is a liquid and the other is a vapour surrounding the 
liquid and solids; the surrounding fluid may as well 
be a second, immiscible liquid. In the following -dis- 
cussion, only configurations in which this fluid is 
simply connected and the solids are bounded by 
convex surfaces will be considered. 

Each equilibrium configuration of the system 
is characterized by the degree of connection, c, of 
the liquid and by the type of the simply connected 
volumes of liquid, or liquid regions (see examples in 
Fig. 1). These liquid regions can be classified in terms 


of two integers n; and n., respectively defined as the 
number of solids completely immersed in the region 
and the number of solids partly wetted or contacted 
by -the region. For example, for an isolated drop 
n=n.=0; for a sessile drop n= 0, n.= 1; for a 
liquid bridge between two solids, n;, = 0, n. = 2; for 
a liquid layer on a solid, n,= 1, n.=0. It is also 
necessary to specify which solids are immersed or 
contacted by the region. Configurations differing only 
in the volume of liquid in each region or in its loca- 
tion, but with the same number and type of regions, 
will be termed similar configurations. In Appendix 
1 it is shown that the number of non-similar liquid 
regions in a system with n solids is 3º, This number 
is reduced to (n + 1) (n+ 2)/2 if the solids are 
treated as undiscernible; for example, a single drop 
on a solid may be counted as n regions or just as 
one region. As the degree of connection, c, is allowed 
to increase, the mumber of configurations rapidly 
increases. For example, for two solids (n = 2) the 
numbers are as follows: c=1: 9 and 6; c=2; 20 
and 13; c=3:34 and 21; c=4: 53 and 32. 

We shall distinguish between two types of equili- 
brium configurations. In normal configurations, the 
equilibrium is stable in all regions and there is no 
contact of interfaces in a region or between different 
regions. As the relative positions of the solids and 
liquid regions are changed, or the volume of liquid 
in a region is altered, starting at a normal configu- 
ration, a critical configuration may occur. This is the 
other type of configurations, examplified in Fig. 2. 
A critical configuration may result when two or more 
interfaces come in contact, such as two drops in 
contact, or when a continuous sequence of stable 


Fig. 2 — Examples of critical configurations which lead to 
contact transitions (a,b) or overshoot transitions (c). The 
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arrows indicate the directjon of displacement of the rods to 
induce the transition 


equilibrium configurations leads to a limiting confi- 
guration which is unstable. A very well-known exam- 
ple is the limiting configuration of a pendent drop, 
just before detachment. Another example, more rele- 
vant to the following discussion, is the configuration 
of a bridge in which a line or point singularity (neck), 
due to contact of the fluid interfaces, appears in the 
bridge when the separation of the two solids reaches 
a critical value. The critical configurations due to 
contact of interfaces can be classified according to the 
interfaces that come in contact. The possibilities are 
LV-LV, LV-SV and LV-SL, where S, L and V denote 
the solid, liquid and second fluid (vapour) phases. 

The critical configurations may be stable or uns- 
table. One of the topics that will be discussed in 
this paper is the stability of critical configurations. 
When a critical configuration is unstable, a new equi- 
librium configuration will result if the system is per- 
turbed. The other main topic is the prediction of this 
final configuration. This is a dificult problem, which 
has not been solved even for the simple case of drop 
detachment, 

Ki is convenient to classify the transitions into two 
classes: contact transitions, when the critical unstable 
configuration is produced by contact of two different 
regions; and overshoot transitions, when the stability 
limit of a single region is reached. 

Examples of phenomena where these transitions are 
relevant are the coallescence and breaking of liquid 
drops on a solid surface and the formation and rup- 
ture of liquid bridges. Related practical areas include 
wetting, drop detachment, adhesion, powder techno- 
logy, emulsion stability and cell fagocitosis. 

This paper contains a detailed analysis of the equi- 
librium configurations and of the transitions at criti- 
cal configurations in a simple system with cylindrical 
symmetry and comprising two equal rods of circular 
cross-section. Gravity will be neglected. The treat- 
ment of this system is greatly simplified by the fact 
that the cross-section of the fluid interface is circular 
ior all types of liquid regions. The analysis of systems 
with more common simple geometries (e.g. two sphe- 
rical particles) is considerably more complex. A simi- 
lar study has already been undertaken for two parallel 
plates with cylindrical symmetry [1]. 


2 — CONFIGURATIONS AND GENERAL EQUA- 
TIONS 


The system includes one or two circular rods of 
radius R and length L, with parallel axes, and a vo- 
lume V of liquid. The system is bounded by two 
parallel plates, at a distance L, such that the angle 
of contact of the liquid with the plates is 90º. All 
liquid regions contact and are bounded by these two 
plates. The 90º angle of contact ensures the required 


cylindrical symmetry of the liquid regions. For exam- 
ple, what will be termed as an isolated drop is in 
fact a cylindrical volume of liquid between the two 
plates, 1.e. a liquid bridge between the plates. This 
artifice in the definition of the system to include the 
two side plates is similar to that previously used in 
studies of pendent drops [2] and of liquid bridges 
between rough plates [3], both with cylindrical sym- 
metry. On the other hand, by reducing the separation 
L we can always avoid problems of instability of the 
liquid regions in relation to configurations without 
cylindrical symmetry [2]. We also note that the con- 
tribution of the plates to the free energy is constant, 
independent of the wetted area in them, which is a 
consequence of the 90º contact angle. 

The contact angle, 8, of the liquid with the rods 
(in contact with the second fluid) is assumed to be 
a characteristic of the system and related to the 
interfacial tensions Ysv, Ys and Y (the fluid interfacial 
tension) by Young's equation 


| tsv — Tae 
cosê = = = E, (1) 


When the quantity E is outside the interval (— 1,1), 
a contact angle cannot be defined and some types of 
liquid regions (e.g. sessile drops) are ruled out. 

In a system with one rod only, the types of equi- 
librium liquid regions are indicated in Fig. la: isola- 
ted drop (1), layer (L) and sessile drop (S). When the 
system comprises two rods, additional types can occur 
(Fig. 1b): liquid bridge (B), rod enclosed in a sessile 
drop (E) and connecting layer (C). The regions L, E 
and C may be symmetrical (if they have maximum 
symmetry) as in Fig. la, b, or non-symmetrical as in 
the examples of Fig. lc. Bridges always admit two 
planes of symmetry. 

The multiply connected configurations, exampli- 
fied in Fig. 2, will be indicated by specifying the 
component regions, for example, [I+S, 21, B+sS, 
E + S. When ambiguity occurs, such as for two S in 
a two-rod system, the notation 2S indicates that the 
two drops are on the same rod, while S+S means that 
they are on different rods. When |E | > 1, the regions 
B, S and E are not possible. 


For each configuration we shall calculate the Hel- 
moltz energy per unit lenght along the cylindrical 
axis of the system. Since we are only interested im 
comparing energies of equilibrium configurations 
with the same volume of liquid and at the same 
temperature, it is enough [4] to calculate the Helm- 
holtz energy terms, A”, due to che interfaces with the 
rods and to the fluid interface: 


A' = Osv tsv + Os TsL += 117 (2) 
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where the QQ are the areas of the interfaces per unit 
length (Q refers to the fluid interface). But Osy + Qsz 
is a constant for a particular system, so that we may 
obtain the Helmholtz energy per unit length from 


A=Y(2-— Eng). (3) 


This equation shows that the value of A for any 
configuration in a one- or two-rod system can be 
obtained by addition of the A values of the compo- 
nent regions. 

In Table I are indicated the main geometrical 
characteristics of each liquid region, including the 
interfacial areas that enter in eg. 3. The definition 
of the geometrical parameters is given in Figs. la, b: 
ris the radius of I, e the thickness of a symmetrical L, 
h the height of S or E, c the radius of C and £ the 
distance between the rod surfaces; p is the radius of 
curvature of the fluid interface in S and B regions; 
& is the wetted arc on a rod, w is the arc subentended 
by the fluid interface in S or B (in bridges, p and w 
have the same sign; » > O means that the centre of 
curvature is in the side of the liquid); finally, n is the 
neck width of a bridge. 

All linear dimensions will be expressed in units 
of R, the rod radius, and denoted by an asterisk, e.g., 


The volume per unit length is expressed in units of 
=Rº?, the volume of the rod per unit length, 


ve=—— (5) 


Finally, the reduced Helmholtz energy, A*, per unit 
length is defined by 


A 


de — ind o 
where 
Q Os 
o— — * oO — E 7 
a 2nR * Us 27R 0) 


are the reduced areas. 

The limits of geometrical possibility of each liquid 
region are also indicated in Table I. The limits for 
L, C and E configurations given in Table I apply to 
the symmetrical configurations. The maximum ?* for 
a bridge with given V* corresponds to n* = 0 for 
0 < 90º and to q =0 for > 90º, the minimum ** 
may be zero or occur for p = 180º when € < 90º. At 
constant ?*, the extreme volumes for B are as follows: 
for € < 90º, minimum n = 0, maximum g = 180º; for 


pise qtas fi Es ira (4) 8 > 90º, minimum & =0, maximum V*= o (for 
Rº R' cd & = 270º — 8). For 0 = 90º and constant V*, 2* may 
TABLE | 


GEOMETRICAL PROPERTIES OF 


Region Reduced radius 


Volume Vº 


| Il 
isolated drop 


L 
layer 


E*(c* 42) 


TVE = p*Hu + 


S 
sessile drop 


+p*sind—q 


C 
connecting 
layer 


E 
enclosed rod 


as for 5 


mV a =4p*sinto 


sing +p*2k(to)— 
=k(p); 
k(0)=204— 
sinZa 


LIQUID REGIONS 


Area solid/ 
liquid 

O] 
sz SL 


Height; other 


) Range of existence 
relations 


r*>0 
E*>0 
O<p<m— 
p* 
c*22+ 2 
h$22 
| 2420; n*20; 
| Ore-mwm=71/2; —T<w<tT/2; 


e<m if 0<m/2; 
Im 
dão REC —0;[0>7/2 


n*=sing + 
+p*(I—costw) 


vary between O and co; and for constant **, V* may 
vary between O and co. Fig. 3 gives the extreme values 
of B for three values of 6. For € < 90º there are two 
curves in Fig. 3, respectively for the maximum ** or 
minimum V*. 

In bridges it may also be required to know the 
force F (per unit length) that has to be applied to 
the rods to keep equilibrium. This force is the sum 
of a surface tension force and a pressure difference 
force [5]: 

sing 


F =2y (cosw — =. » (8) 


It can also be obtained from the slope of the curve 
A*(º*) at constant V* [5]: 


dA* F 
p=" =>— (9) 
de* 27% 
|* 
12 30, | B=0º 
g| “8=150º 
6 
B=0º 
31 
B = 30º 
0 3 6 9 q 5 vº 


Fig. 3 — Limits of geometrical possibility of liquid bridges 
for various values of the contact angle 


3— HELMHOLTZ ENERGY OF CONFIGURA- 
TIONS 


Curves of the Helmholtz energy, A*, as a function 
of V* can easily be prepared for each type of region, 
and for selected values of € (and ** in the case of 
bridges). Examples of such curves for € = 30º, 90º 
and 150º are given in Fig. 4. The curves are -of posi- 
tive slope and concave downwards except for bridges 
with 0 < 90º, in which case a minimum (F = 0) 
followed by an inflection point may occur, as in 
Fig. 4a, 

The value of A* for any multiply connected con- 
figurstion can then be calculated from these curves. 


8 =30" ] 


Bl =0 PEL 


= 90" 


"RR Lv ÕÃr% By 
(e) 
Fig. 4— Variation of MHelmholtz emergy with volume of 
liquid regions for the contact angles indicated. An arrrow 
indicates that a curve continues. A dot indicates the extreme 
value for the region. In (a) the curves E, C are nearly 
coincident; the same occurs in (c) for the curves E, L; im 
both cases E is slightly below 
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Fig. 5 shows examples of the variation of A* in 
doubly connected configurations of fixed total volume 
V*, as a function of the volume in one of the regions, 
for 0 = 30º. These curves can be obtained by a sim- 
ple graphical construction from the A*(V*) curves for 
the component regions, as shown in Fig. 5a. For the 
pars 855 (07:28), I+FI, 2L IFC. I+E and 
S + E the curves are similar in that they have a 
maximum A* for a particular composition of the 
system (i.e. a particular distribution of the total volu- 
me between the two regions); for I + L a maximum 
occurs if V*> 1. For the L region, the value of 
A* for V*t= 0 («= 0) is taken as the one of a rod 
with a film of liquid; the value of A* for the dry rod 
is larger by (cosê — 1). The energy curves for B+I 
and B + S are also of this type for € > 90º. However, 
for 0 < 90º in addition to the maximum there may 
be a minimum A* as in the example of Fig. 5c. 

An important property of the stationary points of 
A* for doubly connected regions is that they corres- 
pond to equal radii of curvature of the fluid inter- 
faces in the two regions. This is obvious from Fig. 5a 
when the two regions are of the same type (21, 2L); 
the proof for the other cases is given in Appendix 2. 
This property is illustrated in Fig. 5c for the confi- 
gurations B+I and B+S (with 2* = 0). 

For the purpose of discussing transitions at critical 
configurations it is also useful to construct curves 
A* (2*) at constant V* and €, such as those in Fig. 6. 
The volumes in each region in multiply connected 
configurations are assumed to be equal. Other confi- 
gurations could of course be included in Fig. 6. The 
curves show the range of * on the assumption that 


Fig. 5 — Helmholtz emergy in doubly connected configura- 
tions as a function of the volumes in each region (8=30º). The 
dashed curves in (a) are fo: the pure regions indicated. In 
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all configurations are symmetrical. The critical con- 
figurations are then reached at the extreme values of 
t*. For non-symmetrical configurations the limits 
of ** are of course different: the critical configura- 
tions may be reached with ** = 0 in I,S+S,L,2L, 
L+I, L+S, C and E configurations (see Fig. 2a). 

From the numerical results examplified in Figs. 4- 
-6 a number of general conclusions can be drawn. 
A condensed notation is used to compare the energies: 
for example nI > I means that n isolated drops have 
higher A* than a single drop with the same total 
volume. 


One-rod configurations: 


ni>I 

nS>sS 

[= 

L>S 

For obtuse 8: L > I. For acute 0, L > I at small 
volumes, but I > L at large volumes. The critical 
volume at which I=L is given as a function 
of 0 in Fig. 7a, b. 


q pa e 


Two-rod configurations: 


6 C>E(C=Eforê->õo0) 

L>E(L=E for ê%-s 180º) 

8. For each of the pairs (C, L), (€, 1, (€, S), (E,S), 
(E, 1) the relative energies depend on V* and 0 
in a way similar to that described above for (L, 1). 
Fig. 7c refers to C=. 


= 


* 


9. The bridges are attractive > 0) at large 


ro 


(c) are also shown the radii of curvature of the component 
regions; the p* curve for B has two branches with a common 
asymptote. The B configuratijons is for ?*=0 


Fig. 6 — Helmholtz energy a function of rod separation 

for fixed volume and contact angle. The dot indicates the 

extreme value of £* for each configuration assumed symme- 

trical. In multiply connected configurations, the volume is 
equally distributed in the component regions 
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separations and repulsive at small separations, 
except for small volumes and 8 < 90º in which 
case they are always attractive (see Fig. 6a). 

10. Among all configurations compatible with a 
given position of the solids and a given volume 
of liquid, the one with smaller A* is always 
simply connected. In particular, the smaller A* 
for configurations with two regions occurs when 
one of the regions disappears. This is true even 
if one of the regions is B (8 < 90º) as in Fig. 5c. 


4 —STABILITY AND TRANSITIONS BETWEEN 
EQUILIBRIUM CONFIGURATIONS 


It was implicit in the previous analysis that the 
equilibrium liquid regions are stable in their interval 
of geometrical possibility, provided they do not come 
in contact with other regions. This is true even if 
there are other configurations with smaller A*, 


The two classes of critical configurations are illus- 
trated in Fig. 2. Fig. 2a, b shows examples of critical 
configurations due to contact of interfaces: in (a) the 
contact is LV — LV and in (b) the contact is SL — 
— LV. Fig. 2c shows examples of critical configura- 
tions that lead to overshoot transitions; these may 
only occur in bridges, in the system under conside- 
ration. 

It is now necessary to find out whether the 
critical configurations are stable or unstable and, in 
the latter case, predict the new equilibrium state that 
results from a transition induced by a small pertur- 
bation of the system. The critical configurations in 
liquid bridges are certainly unstable, simply because 
there is no admissible neighbouring configuration for 
a separation larger than the critical one. In all other 
cases, we shall use a very simple argument to decide 
on stability. The critical configuration is stable if 
there is no other equilibrium configuration with 
smaller free- energy, A*. If such a configuration 
exists, it is not possible to decide whether equilibrium 
is stable or unstable, unless the configuration with 
smaller A* can be considered as resulting from a 
small perturbation of the critical configuration, in 
which case we may conclude for instability. 

In order to predict the final state of a transition 
from an unstable critical configuration, it may be 
necessary to define precisely the type of perturbation 
that destroys equilibrium. For example, in overshoot 
transitions in liquid bridges it is necessary to specify 
the relative displacement of the rods that perturbs 
the system. We shall consider symmetrical perturba- 
tions, in which the two rods are symmetrically dis- 
placed. 

In general, the prediction of the final state in a 
transition will be made by considering the values of 
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A* of the candidate final configurations and invo- 
king what could be termed a continuity argument, 
i.e. the existence of a continuous path for the system 
between the initial critical configuration and the final 
state. Finally, a dynamic argument will also be used 
for contact transitions, related to the pressure in the 
liquid regions in contact. The pressure in the surroun- 
ding fluid is uniform (P,), but in each liquid region 
the pressure, P, will depend on the radius of curva- 
ture, p, Of the interface, according to the equation 
of Laplace; 


= pa, (10) 


p can be positive or negative (in bridges). When two 
regions contact, it is expected that, in the transition, 
the liquid flows from the high pressure to the low 
pressure region, at least initially. The following 
examples strongly suggest that the transition occurs 
in this direction; in some cases the final configuration 
contains regions with the same liquid pressure, but in 
general the high pressure region disappears comple- 
tely. 

We shall now apply these principles and argu- 
ments to assess the stability of critical configurations 
and to predict the transitions from unstable configu- 
rations. We consider in turn systems with one and 
two rods. 


a) One-rod systems 
Only contact transitions occur in such systems. 


l. When two or more I contact, a single I results 
(Fig. 5a). The decrease in A* would also occur 
if the final configuration contained two drops, 
one larger and the other smaller than the initial 
drops. However, it is expected that only one 
drop results, its position being completely deter- 
mined from the invariance of the position of the 
mass centre. It is as if the system would follow 
a continuous path of equilibrium configurations 
with decreasing A* and its final configuration 
corresponds to the minimum A* in the sequence 
of equilibrium configurations, which, of course, 
is not the actual path of the system. The actual 
path could only be determined from a dynamic 
analysis of the transition, and this type of ana- 
lysis is extremely difficult. 

2. When two or more S contact, a S results. Note 
that L cannot form, since the wetted arc decre- 
ases by coallescence of S drops. The final posi- 
tion of the rod and S cannot be determined from 
momentum conservation. 

3. Contact L, I: the final configuration is L or I 


VR TT me e E 


depending on the relative values of the radii of 
curvature in the initial configuration. It is as if 
the energy barrier in the curves of Fig. 5b is avoi- 
ded and the system followed a path of decreasing 
A*, For |E| < 1, when I results, the rod becomes 
dry. However, for E> 1 it is expected that a 
liquid film remains on the rod, because in this 
case the dry rod would give a total energy larger 
than the initial one. 

. Contact I, S: analogous to 3. The final confi- 
guration is I or S as if the system would follow 
a path of decreasing A* (Fig. 5a). If the contact 
is symmetric, the final position of the rod is 
determined. 

- Contact S, S: the final state is a single S for the 
reasons explained in 1. 

. An I drop of reduced radius r* in contact with 
the dry rod can, in principle, produce L or S, or 
be stable. S can only form if [E|< 1 and always 
leads to a reduction in energy. For the forma- 
tion of L, the reduction in energy only occurs 
for combinations r*, E to the left of the 
curve in Fig. 7a,b. Intermediate equilibrium 
states S+I and L+I willin most cases have 
larger energy than the initial I (Fig: 5a, b). The 
argument given above in 3. then suggests that the 
isolated drop is stable in such cases. The excep- 
tion occurs in the formation of L 1f the solid 1s 
sufficiently wettable, It is easy to show that for 
E> (21% + 2)á — r*, which implies E > 1, the 
initial energy (1 + dry rod) is larger than the 
maximum energy of the intermediate L + I con- 
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(a) 
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figurations. We expect therefore that in such 
cases the drop will be unstable leading to the 
formation of L; otherwise I will be stable. The 
engulfment of particles by a liquid was discussed 
by Neumann [6], but he considered only the 
condition of energy decrease. The present ana- 
lysis suggests that this necessary condition is not 
in general a sufficient one. 


- In all cases the degree of connection decreases 


in the transition; the higher pressure region 
disappears. 

When the radii of curvature of the two regions 
in contact is the same, the configuration is also 
unstable but the final configuration depends on 
the way the system is perturbed. For example, 
in the contact of I and S with equal radii, the 
following perturbation can be envisaged. A small 
amount of liquid leaves I, for example; then the 
radius of I will become smaller than that of S, 
and more liquid leaves I. And conservely, if 
the volume of I increases slightly. In both cases 
the configuration is unstable. 


Two-rod systems 


We consider separately contact transitions and 
overshoot transitions. 


Contact transitions 


In addition to the conclusions drawn for one-rod 
systems we have: 


ve 


AS SS SS SS SS Sã SS + 


0 % 60 90º 


Fig. 7 — Critical volume, V*, for which the energy of L and 
I (in a and b) and of € and I (in c) are equal, The more 
stable region in indicated in each case 
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10. 


ll. 


12. 


13. 


14. 


15. 


L6. 


17. 


18. 
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Contacts 1, C and I, E: the final configuration 
depends on the relative pressures in the initial 
state. Point 7. applies with the restrictions dis- 
cussed in 3. as regards the case of E >. 


Two sessile drops symmetrically placed on each 
rod (S + S) contact to give B (or C if the volume 
overshoots the maximum for B). See Fig. 6. 
Contacts B, I and B, S. When the A*(V*) curve 
has no minimum (e.g. Fig. 5b), the final confi- 
guration is simply connected and can be pre- 
dicted as discussed in 3. and 6. If the A*(V*) 
curve shows a minimum (acute 8, not too small 
volumes) as in Fig. 5c, the final configuration 
is probably the one of smaller energy that can 
be reached with decreasing A*. This means that 
for some initial configurations, a doubly connec- 
ted configuration will result in the transition, 
with a composition corresponding to the mini- 
mum. The special configuration at the minimum 
is stable. This is confirmed by an argument 
such as the one in 8.: the radii change in such 
a way as to re-establish the original volumes. It 
is interesting to note that there is no other con- 
figuration (with other regions) of smaller A*. 
When S is contacted by the second rod, the 
transition is probably to B for symmetrical con- 
figurations (Fig. 6). For non-symmetrical confi- 
gurations it may happen that B cannot result, 
such as in Fig. 6d-f for ** > 0, In these exam- 
ples, there is no alternative configuration with 
lower A*; the contacted S is probably stable. 
When I is contacted simultaneously by the two 
rods, similar conclusions apply. B results in the 
symmetrical case, but otherwise B may have 
higher energy. S+S or 2L are ruled out as 
final configurations; they should have equal 
volumes in the two regions, in which case A” 
would increase (Fig. 6). 

When L contacts the second rod, similar con- 
clusions to 12. and 13. apply. If L is not centred 
and the contact is for small 2*, a B may be 
forbidden, as in Fig. 6d. 

When two L contact, similar conclusions apply. 
B results under symmetrical conditions (equal 
volumes in 2L and centred L) but may be for- 
bidden under asymmetric conditions and diffe- 
rent volumes in each L. Alternative possibilities 
are Cor L. 

When L and S contact, the result can be B, L 
or S, depending on the volumes of L and S. 
When two L (or S+ S) contact simultaneously 
I, under symmetrical conditions, the result is B 
(Fig. 6). 

In non-symmetric contacts, C may form instead 
of B if the total volume is too large (acute 8). 


For example, if I is very large and contacts 
simultancously the two rods at a short distance, 
C will result (Fig. 6b). 

19. Some contacts may be geometrically forbidden. 
This is the case for the contact of a large 1 
with a bridge of negative curvature. These for- 
bidden contacts do not alter the conclusions in 
l1. 

20. If the maximum separation for C is reached, 
leading to contact of the rod with the surface of 
C, B results at least under symmetrical conditions 
(e.g. Fig. 6b). Otherwise, B may be forbidden 
and E probably occurs (Fig. 6b). It is however 
likely that this critical configuration is stable 
(see below). 

21. If the maximum separation for E is reached 
symmetrically, B results (e.g. Fig. 6d). In non- 
-symmetrical situations B may be forbidden and 
it is difficult to guess the final configuration. 
For example, in Fig. 6b there is no configuration 
with smaller A* than E as 4* —» 0. This suggests 
that the E configuration is stable, and the 
immersed rod can only leave the liquid if a 
finite amount of energy is supplied. 


b2) Overshoot transitions 


22. If the maximum separation for B is overshoot 
symmetrically, the final configuration is I for 
O > 90º (Fig. 6e, ) or S +S with equal volumes 
for é < 90º (Fig. 6a). A non-symmetrical pertur- 
bation may lead to other final configurations. 
23. If the minimum separation for B is overshoot 
(0 < 90º, large volumes) a C results (Fig. 6b). 


It is concluded, therefore, that liquid-liquid con- 
tacts are unstable but other contacts are stable if a 
contact angle can be defined (| E | < 1). 

When a single region results from a transition, 
the invariance of the mass centre may be sufficient 
to determine the position of the region, which is 
necessarily at rest. In such cases, it is as if the region 
that survives were attracted by the other region. 
If two regions result from a contact, their positions 
and individual momenta remain undetermined. 

As referred above, when a B is formed in a tran- 
sition, forces must be applied to the rods to keep 
equilibrium. If such forces are not applied, the final 
configuration can be quite different, in general with 
a change in *. For example, when S+S contact 
and a B is formed, if forces are not applied to the 
rods thes: will move in a direction that depends on 
whether the force is attractive or repulsive. In the 
first case the rods will eventually contact, but if the 
force is repulsive the rod separation may increase 
initially (the force is always attractive at large dis- 
tances). 
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5 — CONCLUDING REMARKS 


The method developed in this paper to predict 
the transitions at critical configurations avoids the 
description of the irreversible transformation that 
occurs and is essentially based on a comparison of 
the properties of the initial and final configurations. 
The merit of the method is that it rules out some 
final configurations and leads to simple predictive 
rules in the case of contact transitions: the final 
configuration is the one of smaller A* that can be 
reached from the initial configuration following a 
path of decreasing A* in the diagrams of the type 
shown in Fig. 5. Finally, the analysis shows that 
there are stable critical configurations which can 
only be destroyed if mechanical energy is supplied 
to the system; these are the configurations resulting 
from the contact of a dry solid with a drop, or of a 
wetted solid with a liquid surface, when |[E|</ 1. 

As pointed out in the Introduction, when a critical 
configuration is reached in the system discussed in 
this paper, a line singularity occurs in the surfaces. 
When the configuration is unstable, a transition 
takes place to a new configuration with no line singu- 
larities. Energy terms associated with these lines 
should then be included in the free energy of the 
system and, for the sake of coherency, the energy 
contributions of all lines of contact between two 
different surfaces (i. e., the solid-liquid-second fluid 
lines of contact) should then also be included [e.g. 7]. 
This would make the analysis much more compli- 
cated and general conclusions on stability and transi- 
tions could hardly be drawn. However, it is expected 
that the relative contribution of the line energies 
will be negligeable when the linear dimensions of the 
system are not too smal. The conclusions drawn above 
on the stability of critical configurations and on the 
final state in a transition should therefore be valid 
in such systems, but not in those with small dimen- 
sons when the contribution of line energies is not 
negligeable. 


Appendix 1 — NUMBER OF CONFIGURATIONS 


The number of simply connected liquid regions 
that may occur in a system with n solids (convex 
surfaces) is obtained as follows. Each region is charac- 
terized by the integers n;, n. (> 0): the first indicates 
the number of solids immersed in the region and the 
second the number of solids partly wetted by the 
region. The maximum value of k = n; + n. is n, and 
all combinations (n;, n.) such that n;+n. <n can 
occur. For example, for n = 2, the pairs (n;,n,) are: 
(0,0) (1,0) (0,1) (2,0) (1,1) (0,2). For a particular pair 
(n;, nc.) the number of distinct regions is k + 1 if the 
solids are treated as undiscernible and (k)- (2) other- 


wise; (7), for example, denotes the number of selec- 
ting k objects out of n. For example, for n,=3, 
ne =1, there are ()=4 distinct regions involving 
any selected 4 solids. If the number of solids is n, 
the 4 solids can be chosen in (1) ways. The total 
number of regions in the first case is therefore 


(n + 1) (n + 2) 


N=1+2 et (nr = 3 


In the second case, the total number is 


n k n 
N=5 [3 O | =5 200) =3 
0. j=0 k=0 


k= 


Equations for the number of configurations with 
higher degrees of connection, c, are hard to obtain. 
The values given in the text for c=2 were found 
by writting down all the possible configurations. 


Appendix 2— STATIONARY POINTS IN A*(V*) 
CURVES FOR DOUBLY CONNECTED CONFI- 
GURATIONS 


Consider first the configuration S + I. The total 
volume (cf. Table 1) 


Vº=r"+ — [e* (e*w + sin 0) — &] (A.2.1) 


is fixed. r* and p* are the radii of curvature of I and 
S, respectively: 


in 
p'P=-—— ;vu=0+6: (A2.2) 


The total A* is 


A*=r' + — (up — cos 8) - (A2.3) 
From dV* = 0 we obtain 
2mrº dr + Ea +sin0 + 
E 2) 
ia [o = (A2.4) 
sin 8 | 
If A* is stationary, dA'=0, and 
a) | 
ndr* + E + e (o* — cos o | de' = 0, (A2.5) 
| sin 8 | 
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Combining A2.4-5 gives 


sinão 
Ze*'w + sinb + —— (p*2— 1) 
sin ê 
2r'=—D—————— -——————— (A2.6) 
sin? (» 
q + —— (p* — cos 8) 
sin & 


Simple calculations, using eqgs. A2.2, give the final 
result 


r =p". (A2.7) 
The stationary point corresponds to a maximum, as 
explained in the text. 

Consider now the configuration B+I. The equa- 
tions are 


l 
dd ade dieeé is sin w sin & — 2(6 — 


— sin &cos &) + 2e*' (w — sinwcosw)] (A2.8) 
ma 

A*=r* + EE (p*w — & cos 0) (A2.9) 

with 
0—-w+ & =7/2:; do = dd (42.10) 

and 

"M : | 

=p'sinw+cosb-—l (A2.11) 
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so that 


(— sin & + p* cosw) do + sinw p* = (A2.12) 
The conditions dV*=0 and dA*=0 combined with 
egs. A2.10-12 give 


=p". (A2.13) 


The stationary point may be a maximum or a mi- 
nimum. 

Combining the previous results, it is easy to con- 
clude that for B+S configurations the stationary A” 
corresponds to pj=p( - 

The same property holds in all other cases, such 
as 2I, 25, 2L, L+I, L+S, I+E, I4C and a formal 
proof can be made on the same grounds. 
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ANÁLISE DOS METODOS DE CÁLCULO 
DE CARGAS TÉRMICAS ATRAVES DE 


ESTRUTURAS 
L.F. RORIZ* 


RESUMO 


Neste artigo são discutidos os principais métodos de 
determinação de cargas térmicas através da estrutura de 
edifícios, sendo analisados os erros teóricos inerentes a 
esses métodos. 


SUMMARY 


An analysis of the three main methods for the 
evaluation of the thermal load through structures is 
presented. The methods recommended by Carrier, IHVE 
and ASRAE respectively, use one of the two fictitious 
temperatures: the sol-air and the equivalent temperatures. 
The errors occuring using each method are pointed out, 
and the general conditions where each method can be 
used as well as those where each method is advisable 
are also given. 

Carrier method has only the advantage of beeing 
speedy but its use is not advisable as it can lead to large 
errors in the calculated load. The use of the ASHARE 
method is advisable but only for the cases presented in 
the tables. From all methods the most general with 
good precision is the IHVE method. 


| - INTRODUÇÃO 


Os métodos discutidos neste artigo são os de utiliza- 
ção mais generalizada Estes métodos descritos em [1], 
[2] e [3] utilizam uma das temperaturas fictícias usual- 
mente designadas por temperatura sol-ar (T,) e por 
temperatura equivalente (T. ). 

A temperatura equivalente é a temperatura fictícia a 
que o ar exterior se deverá encontrar de tal forma que o 
fluxo de calor que atravessa a superfície interior da 
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Original recebido para publicação em 12/09/85. 


estrutura, ie, que entra no espaço a climatizar devido à 
radiação e à diferença real de temperaturas entre o ar 
exterior e interior é igual ao que entraria nesse espaço a 
climatizar, em regime estacionário, não havendo radia- 
ção e encontrando-se o ar exterior a essa temperatura 
fictícia. A temperatura sol-ar é a temperatura fictícia 
a que o ar exterior se deverá encontrar de tal forma que 
o fluxo de calor que atravessa a superfície exterior da 
estrutura devido à radiação e à diferença real de tempe- 
raturas entre o exterior e o interior, é igual ao que atra- 
vessaria essa mesma superficie, não havendo radiação e 
estando o ar exterior a essa temperatura. 

Embora a carga térmica através da estrutura de edi- 
fícios seja do tipo sensível e do tipo latente, todos os 
métodos que são analisados neste artigo apenas consi- 
deram a carga sensível. Todos os métodos desprezam 
também o efeito dos bordos, considerando as superfi- 
cies estruturais independentemente. Esta forma de pro- 
ceder é aceitável pois o erro introduzido é desprezável. 
Os erros que serão apontados nos diversos métodos 
referem-se apenas ao método em si, não sendo em qual- 
quer dos casos considerados os efeitos da transmissão 
de massa quer no valor de carga latente quer no valor 
da condutividade térmica aparente a utilizar. 

Dos três métodos que serão analisados, os que se 
baseiam, na temperatura equivalente [1] e [3] são de 
utilização mais simples e rápida que o que se baseia no 
conceito da temperatura sol-ar [2]. Para os distinguir, 
os métodos serão designados por métodos Carrier, mé- 
todo IHVE e método ASHRAE. 


2 - MÉTODO CARRIER 


Neste método a carga térmica é calculada através de 
Q=-AK, AT, (1) 
onde K, é o coeficiente global de transmissão de calor 
da estrutura considerada e A é a área da estrutura. 
Os erros relativos a este método foram anteriormente 
discutidos pelo autor [4]. 
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Como então foi referido o fluxo real de calor que 
atravessa a estrutura durante 24 horas terá que ser igual 
ao que atravessaria a estrutura em regime estacionário 
não havendo radiação e estando o ar exterior à tempe- 
ratura T.. 


uh ps 
f mê ro Res Aa + E 15 = Ay pm ++ Tim = Ea (2) 
sendo AT... a diferença de temperatura equivalente à 


eg; E 
hora i para a parede de espessura j. Este valor da tem- 


peratura fictícia é independente da espessura da parede 
e apenas função da orientação e inclinação da superfi- 
cie em estudo. No entanto tal não se verificar nas tabe- 
las apresentadas em [1]. As diferenças entre os valores 
máximos e mínimos de A [., são apresentadas na Tabe- 
la 1, já que tais valores não foram apresentados em [4]. 
Os valores apresentados foram deduzidos a partir das 
tabelas originais emºF. Da observação da última coluna 
da tabela é fácil deduzir que apenas os valores de AT, 
relativos a terraços em sombra podem ser considerados 
correctos. Para paredes e terraços ao sol o erro no 
valor de AT. é de t+0.6ºC e de +2.2º%, respectiva- 
mente. Por razões de simetria em relação às 12 horas 
solares, valores de T., para paredes simétricas em rela- 
ção à direcção sul deveriam ser os mesmos. Embora 
para o valor médio (AT.,) este erro seja pequeno 
(+0.1ºC) como se pode deduzir da tabela, a utilização 
dos princípios de simetria e da independência da espes- 
sura da parede quando aplicadas simultâneamente a 
paredes E-W resulta numa diferença máxima de AT, 
de 2.67ºF, sendo portanto este o erro real neste tipo de 
paredes. Para as restantes orientações o erro obtido é 
igual ao observado na tabela 1, pelo que os erros reais 
são para paredes NE e NW de 1.71ºF e para as paredes 
SE e SW de 2.13ºF. 


Tipo de a À AA ecitmas) E DE onigenin) 
estrutura " f z E) 


Parede 


Terraço 
ao Sol 


Terraço 
à sombra 


TABELA 1 


Em termos percentuais o erro nos valores tabelados 
de AT., pode ser superior a 100%. 

A utilização das fórmulas de correcção propostas 
em [1] e os erros introduzidos nas diversas correcções 
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foram também discutidos em [4]. 

Estes erros situam-se duma forma geral dentro do 
erro com que são apresentadas as tabelas excepto para 
estruturas ligeiras (telhados, paredes finas, etc.). Neste 
tipo de estruturas e principalmente em telhados, os erros 
inerentes a este método podem ser bastante elevados. 
Nestas condições a utilização deste método não é pois 
aconselhável. 


3 - MÉTODO IHVE 


Neste método a carga térmica à hora H + (9 é dada 
por 


Quis = KA [Ta To *t (Tua = T)] (3) 


sendo T,, a temperatura sol-ar média diária, 2 o desfa- 
samento ou atrazo e f o factor decremental ou de amorte- 
cimento. 

Esta equação é deduzida a partir da equação geral 


O=KAIT-T: 'z reco (i-Wo-m)] (4) 
e . 

em que Ta - Tu - X, cocs(nmt / 12 -'Y,)e em que f, 

e , são o factor decremental e o atrazo referente à 

nésima harmónica. 

A utilização da equação (3) apenas considera a har- 
mónica principal. O erro desta aproximação depende 
portanto da aproximação com que uma única harmó- 
nica consegue descrever a evolução da temperatura da 
superfície exterior ao longo do dia. 

Como o período de oscilação (ou mais própriamente 
dois semi-períodos diferentes) influenciam os valores do 
atrazo e do amortecimento, estes valores deveriam ser 
determinados em função da orientação e inclinação da 
estrutura. No entanto de acordo com [2] o factor decre- 
mental é considerado como função da espessura e den- 
sidade do material, enquanto que o atrazo é conside- 
rado como dependente exclusivamente da espessura do 
material. Esta simplificação introduz um erro bastante 
pequeno no valor do amortecimento, (< 10%), mas tal 
não acontece com o tempo de atrazo. Para paredes NE 
os semi-períodos são cerca de 4 e 20 horas e nestas 
paredes os atrazos calculados de acordo com [2] dife- 
rem em mais de 50% dos atrazos reais produzidos. Para 
paredes SE este erro pode ser da ordem de 30%, como 
se pode ver na Fig. 1. Nesta, os valores do atrazo indi- 
cados referem-se a paredes de cor escura, para as prin- 
cipais orientações. Foram obtidos através da solução 
numérica da equação do calor considerando uma varia- 
ção térmica diária de 24 a 35ºC, a radiação solar média 
do dia 21 de Junho e a velocidade do vento de 12 km/h. 
A temperatura interior considerada foi de 27ºC e as 
propriedades das paredes eram: p = 2600 Kg/m', k - 
0.72 Kcal/h mºC, C, = 0.20 Kcal/KgºC. As espessuras 
estudadas foram de 0.105, 0.180, 0.285 e 0.345m. O 


erro total na utilização do método, é a soma dos dois 
erros anteriormente citados. Consoante a orientação, estes 
podem tender a acumular-se ou não. Como se depreende 
da figura o erro final total é mínimo para paredes SW, 
We S (e consequentemente para terraços) e máximo 
para paredes NE. No entanto este erro tende a aumen- 
tar com a espessura da parede. No caso da parede de 
0.345m, onde o valor do atrazo calculado difere, 
para diversos casos, em mais de uma hora do atrazo 
real, o amortecimento (ATim/ AText) é inferior a 4%. 
Para o caso de paredes de 0.285m de espessura o amor- 
tecimento é de 6 a 7%. Nestes casos um erro de 2 ou 3 
horas no valor calculado do atrazo, não afectará signi- 
ficativamente a carga térmica através da estrutura 
calculada. 

Existe ainda outro erro adicional no cálculo do valor 
do atrazo. Este erro não resulta do método própria- 
mente dito mas sim do facto de não se ter acesso à hora 
das temperaturas máxima e mínima da superfície exte- 
rior da parede mas apenas à variação da temperatura 
sol-ar. Há pois que considerar ainda o atrazo devido à 
convecção. No entanto para as diferentes orientações 
este valor varia apenas entre 0.3 e 1.2 horas ou seja, 
para efeitos de cálculos práticos de 0.5 ou 1 hora, con- 
soante o semi-período da oscilação é inferior ou supe- 
rior a 12 horas. 


a. 
semi -perio do da 


oscilação (A) 


4 - MÉTODO ASHRAE 


Este método calcula a temperatura sol-ar de diver- 
sos tipos de estruturas, a partir da qual e utilizando 
funções transferência, calcula o fluxo de calor por uni- 
dade de área que penetra no espaço a climatizar. Com 
base no valor do coeficiente global de transmissão de 
calor da estrutura são determinados os valores da dife- 
rença de temperatura equivalente, designados em [3] 
por diferença de temperatura de carga. 

Desta forma a carga através da estrutura é neste 
método obtido utilizando a equação [1]. Os valores 
foram calculados para uma variação térmica diária: 74 
a 95ºF (23.3 a 35ºC), temperatura interior de 78 ºF 
(25.6ºC), latitude de 40º, mês de Julho, considerando as 
superfícies de cor escura. 

Lógicamente a maioria dos erros inerentes aos 2 
métodos anteriores são evitados. Uma análise análoga à 
referida na discussão do método Carrier é apresentada 
na Tabela 2. Desta forma e utilizando as relações de 
simetria de orientação em relação ao sul o erro máximo 
que é permissível deduzir é <0.2ºF. 

No entanto alguns erros persistem devido às hipóte- 
ses formuladas no modelo utilizado para a variação das 
condições exteriores. De todos o mais importante foi 
considerar que os valores máximo e mínimo da tempe- 


s e Figura | 
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X valores correctos considerado o atrazo como fun- 
ção da orientação, mas incorrectos utilizando o mé- 
y todo IHVE 
: S valores correctos usando o método IHVE, mas in- 
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= orientação 
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pe 
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ratura exterior ocorrem sempre às 15.00 e 05.00 horas 
respectivamente. Dado tal não acontecer, o erro intro- 
duzido no amortecimento e sobretudo no atrazo da onda 
térmica variará com o mês. Os erros por se ter despre- 
zado o efeito do vento e o facto de se considerar radia- 
ção difusa apenas entre o nascer e o pôr do sol quando 
na realidade ela existe entre os dois crepúsculos astro- 
nómicos, são pequenos. O erro cometido no 1.º caso é 
para a grande maioria das estruturas inferior a 5%. O 
erro cometido no 2.º caso apenas reduz ligeiramente os 
valores de cargas térmicas que correspondem às horas 
entre os crepúsculos e o sol atingir o horizonte. A gran- 
de desvantagem deste método é que não é geral e está 
limitado aos casos que foram estudados. No caso de ser 
necessário calcular a carga térmica em paredes e telha- 
dos diferentes dos tabelados não é possível um cálculo 
correcto da carga térmica, mas apenas obter uma esti- 
mativa do valor. 

Quando as condições de projecto forem diferentes 
das que foram consideradas para obtenção dos valores 
de AT. tabelados, há erros que são introduzidos utili- 
zando as fórmulas de correcção apresentadas em [3]. 
Destas, a correcção à variação da temperatura exterior, 
à temperatura interior e à radiação solar durante o ano 
são análogas às do método Carrier que foram já discu- 
tidas em [4]. Para a correcção à cor é apenas apresen- 
tado o valor de 0.5 para cor clara o que é demasiado 
“geral” para todo o tipo de cores. Além disso esta cor- 
recção afecta o valor da transferência de calor apenas 
devido à diferença de temperaturas que é independente 
da cor. Desta forma a utilização de diversos valores 
para o factor de correcção da cor não eliminaria total- 
mente o erro introduzido nesta correcção. 


Orientação | Tipo de Mt | HE..-mM 
| | > im ii 


NH 


E 


z 
E 
on> 


TABELA 2 
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5 - CONCLUSÃO 


Pelo que anteriormente foi dito, dos métodos refe- 
ridos, o do IHVE é aquele que sendo totalmente geral 
oferece maior precisão no valor da carga térmica calcu- 
lada. Para casos gerais a sua utilização é a mais aconse- 
lhável. Embora o método Carrier possa ser utilizado 
se as estruturas não forem ligeiras ou quando a carga 
térmica através da estrutura tiver um peso bastante pe- 
queno na carga térmica total, é aconselhável o uso do 
método IHVE dado existirem tabelas de I, para pare- 
des, terraços e telhados (sendo considerados diversos 
tipos de cobertura) que tornam este último método bas- 
tante expedito. No entanto antes da sua utilização con- 
virá ao projectista verificar se a estrutura que está a 
considerar está tabelada em [3]. Neste caso a utilização 
do método ASHRAE é aconselhável. 
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ABSTRACT 


Á very simple approximate formula for the current 
in a horizontal wire above ground is derived by using 
the quasi- TEM approximation. This formula can be 
interpreted as the current in a transmission line equiva- 
lent to the structure. 


RESUMO 


Utilizando a aproximação quase-TEM deduz-se uma 
fórmula aproximada muito simples para a corrente num 
fio horizontal em presença da terra. Esta fórmula pode 
ser interpretada como a corrente numa linha de trans- 
missão equivalente à estrutura. 


Consider an infinitely long thin wire of radius a, 
parallel to the z-axis located at height b above ground 
and excited by a delta-function generator defined by 
Vô(z) exp(jwt). Let n be the ground refractive index given 
by n=(e,/6,50,/(we,))” with -m/4 < arg (n)<0, where 
e, and o, are the ground dielectric constant and conduct- 
Ivity, and €, is the dielectric constant in free space. Let 
k, = w(u,6,)“ and Z, = (u,/ 65)? be, respectively, the wave 
number and the impedance in free space. The current 
distribution in the wire is given by the integral [1] 


as 2k,V [ N(ah) eJÃZ 
DZ )p DA) 


where T is the real axis, N(ah) = HS(ah) . 7 ah/2; and 
D(h) = hº [T(h) + Q(h)] + Kº [P(h) — Q(h)] (2) 


dA (1) 


Original recebido para publicação em 23/9/85. 


where 


-12bu 
Je a ds 


no. eJ2bu, E Ta 
EU) = T).. u, tu, 4 , QU = T).. níu,+u, | 


T(h) = H,º (ah) - H,º (2bh) 


u =(h-s)2,u,=([(n'-Dki+tu,) with Im(u,), 
Im(u,) <0 for real À and s, h is the transverse wave 
number given by h = (k;- A”) with Im(h) <0 for real À, and 
Hº' is the Hankel function of the second kind and order m. 
The integrand of (1) does not make it possible to derive 
an exact expression for Lin terms of known functions. 
However we shall see that the quasi- TEM approximation 


[2,3] of function D yields an approximation for the current 


in terms of exponentials. 

Neglecting higher modes and choosing in the A-plane a 
contour FT, indicated in Fig. 1.a such that none of the three 
poles corresponding to the three principal modes [2,4,5] 
belongs to the region between contours TT and [j, 
contour Pin (1) may be replaced by contour Ii forz >. 
Furthermore, contour I, (and consequently its image P, 
in the h-plane (Fig. 1.b) has been chosen to satisfy the 
following conditions: 

(1) — The quasi- TEM approximation for D is valid for h 
belonging to Tp. 

(1) — T; and L are steepest-descent paths (SDP) of 
exp(-jAz). 

As it is known [2,3], the quasi-TEM approximation is 
good for |bh| < | except in a region of the complex plane 
corresponding to h close to the branch point [4] h, = , 
(nº + 1)'? - It can therefore be concluded from Fig. 1.b 
that PG is one of the possible contours satisfying condit- 
10n (1). 

Condition (ii) implies that in arcs | and I, where the 
quasi-TEM approximation ceases to be valid (|bh|> 1), 
exp(-jAz) takes values that are sufficiently small for the 
error in the integral of (1) to be neglible when we replace 
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Fig. 1 Integration contour for the quasi-TEM approxima- 
tion. 


L=rt+ERTO N=r+Lh+ E, 


the integrand by its quasi-TEM approximation. This 1s 
possible except if k,z < 1. 

Before using the quasi-TEM approximation of the 
integrand of (1) in order to obtain an expression which, as 
already seen, will approximate the value for current 1, we 
will express function D in the form [6] 


-(hó-K)/Y(h)] (3) 


D(h) = 


where 


z(h)= , [TC + Pe] 


Y(h)=— 


Replacing D in (1) by its quasi- TEM approximation, 
and considering that N (ah) = | for |a h| < 1, we obtain I(z) 
= [o(z) where Ig is given by 


.—. “o exp(JAZ) 
= — => dA 4 
e L. N+ZaYo e 
where 
Zo = Uo)=; + [in(2b 2) +2) So ds 

2Q 

2bs -i 
Po =Y(0)=) pin 

20 
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and uq =V (nº-1)k?-s? . It should be noted that Zq and 
Yo have dimensions respectively of impedance per unit 
length and admittance per unit length [6]. By integrating 
(4) we finally obtain 


Igj=D e fPolz! (5) 


where 


2.= Za and Bo= Va 
are respectively the characteristic impedance and the phase 
constant of an equivalent line of parameters Z, and Yç. 
(5) shows that the quasi-TEM approximation for the 
current is identical with the current which would be obtai- 
ned in the equivalent line referred to above when excited 
by a generator of voltage V. 

In Fig. 2 we give the error associated with the quasi- 
- TEM approximation for the current (5) for the parame- 
ters previously chosen by Carpentier and dos Santos in 
[5]. From the comparison of the results presented in Fig. 2, 


errors, mãÃ 


0.25 
E t=48 MHz 
E 
EM 
0.15 
0.10 
D.05 
1.25 5.0 20.0 80.0 


Fig. 2 Errors for various approximations of the current. 


A=II-L| 4,=II-6| A=[I-L| Aq=II-L] 


| : contribution of the improper mode [2,5] 

I, : contribution of the two proper modes [4,5] 

l. : contribution of the three modes in the form of error 
functions [5] 

ly : the quasi-TEM approximation for the current | 


e/6,=15, 0,=0.01S/m,a=0.0025m,b=Im 


it can be concluded that the quasi-TEM approximation is 
normally more accurate than the one proposed by Chang 
and Olsen [4] (contribution of proper modes only) or than 
the contribution of the improper mode only [5]; but it is 
normally less accurate than the approximation obtained 
by the contribution of the 3 modes in the form of error 
functions [5]. 

It can thus be concluded that, from a practical point of 
view, one can use (5), whose accuracy is surprisingly good 
and which is extremely simple in comparison with repre- 
sentations in terms of modes. 

Finally we would point out that, for very small values 
of k,z, Ly no longer constitutes a good approximation for 
I, owing to the fact that in this case the contribution of 
SDP arcs for |b h| >1 cannot be neglected in the integral 
along T,. On the other hand I and ly have different 
asymptotic behaviours when k,z tends to infinity. This is 
due to the fact that, while the quasi-TEM approximation 
for D is a second degree polynomial in h, D has an asymp- 
totic behaviour in the vicinity of h = o of the type [4] 
C,+C, hºIn h where C, and C, are constants. 
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A TEMPERATURA COMO UM CONCEITO 


TEMPERATURA KELVIN 


RODRIGO DE ABREU 


Instituto Superior Técnico, da UTL e Centro de Electrodinâmica, do INIC 


SUMÁRIO 


O conceito de temperatura Kelvin surge habitual- 
mente sem que se ponha em causa o Primeiro Princípio 
da Termodinâmica. A revisão conceptual do Primeiro 
Princípio feita pelo autor leva a que, dado não ser pos- 
sível atribuir com generalidade sentido físico à expres- 
são de Clausius dE=dW+ dQ, seja necessária a cons- 
trução do conceito de temperatura Kelvin por outra via 
diferente da consagrada por longos anos de tradição. 
Para tal prescindiu-se do Primeiro Principio e do Prin- 
cipio Zero da Termodinâmica. Isto é possivel dada a 
existência duma relação necessária entre energia, entro- 
pia e variáveis de deformação. Uma vez afirmada esta 
necessidade independentemente da noção de tempera- 
tura empírica, adquire-se grande simplicidade e econo- 
mia de conceitos. 


SUMMARY 


The Kelvin concept of temperature generally arises 
without calling in question the First Principle of Ther- 
modynamics. Since it is not possible to attribute a gener- 
alized physical sense to the Clausius expression 
dE =dW+ dQ, the conceptual revision of the First Princ- 
iple carried out by the author leads to the necessity of 
constructing the Kelvin concept of temperature in a 
way different from the one consecrated by long years of 
tradition. To this end, the First Principle and the Zero 
Principle of Thermodynamics have been dispensed with. 
This is possible owing to the existence of a necessary 
relation between energy, entropy and deformation var- 
iables. Once this necessity has been affirmed indepen- 
dently of the notion of empyrical temperature, a great 
simplicity is arrived at, as well as an economy of concepts. 
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INTRODUÇÃO 


Considera-se habitualmente que o conceito de tem- 
peratura é essencial no estabelecimento das relações fun- 
damentais da termodinâmica. 

Mostra-se neste artigo a possibilidade de não se con- 
siderar inicialmente a temperatura (Princípio Zero da 
Termodinâmica) estabelecida que seja a necessidade de 
uma relação entre a energia e a entropia, conforme se 
mostrou em artigo recentemente publicado nesta revista. 

O conceito de temperatura deixa, portanto, de ter 
carácter essencial, passando a ser parâmetro que deriva 
das grandezas energia e entropia. De facto, num sis- 
tema finito, se a energia e a entropia forem constantes 
para diversos valores das variáveis de deformação a 
temperatura também é constante, dado que ao longo da 
intersecção de duas superfícies E=c-É eS=çc tse 
verifica que Ea =€ te (Tese de Doutoramento, I.S.T., 


U.T.L., 1983). Neste artigo mostra-se, no entanto, que é 
possível adquirir com facilidade e com menos abstrac- 
ção o conceito de temperatura, recorrendo-se ao con- 
ceito de “Fonte de Calor” que se define pelos valores 
constantes das densidades de energia e de entropia. 

A demonstração aqui apresentada não sofre de peti- 
tio principii dado que o segundo princípio não se baseia 
na variável entropia introduzida a partir dum conceito 
de temperatura empírica. Este é um dos aspectos essen- 
ciais deste trabalho. 


A TEMPERATURA COMO UM CONCEITO 
DERIVADO DA ENERGIA E DA ENTROPIA. 


Consideremos um sistema de energia E função da 
entropia e do volume. 

O sistema pode ser decomposto em dois sub-sistemas 
satisfazendo as seguintes condições: 


FIG. 


E, = E,(S,, V,) (E,- energia do sub-sistema 1). 

E,=E,(S,)(E,-energia do sub-sistema 2; V;=c—). 

E = E (S,V) (E- energia do sistema constituído 
pela associação dos sub-sistemas | e 2). 


E=E +E,, 

S=5S, + Soy 

V=V,+v, 
e 

Vv,=c 


O sub-sistema 2 é considerado “infinitamente gran- 
de” de tal modo que os acréscimos de Energia por uni- 
dade de volume são nulos e deste modo um sub-siste- 
ma | de 2 (Fig. 1) tem Energia e Entropia constantes 
(E =c Ae S, = + e v, = c!f). Designa-se, habitual- 
mente, o sub-sistema 2 por “Fonte de Calor” ou “Reser- 
vatório de Calor”. 

Consideremos uma variação dS, e dV,, dS,e dV, — 


— temos que: 
dE 
dE, = 350 98 + gy idvo (1) 
e 
0E, | 
dE, = 3, dS» (V;=c Je, (2) 
Como 
E=E, TE (despreza-se um termo (3) 
É,, de interacção) 
tem-se 
dE= dE, + dE,, (4) 
Como 
0É, 
9V, Ma hi P,, (5) 


(se o Bras ae | for o Sistema dado não existir O 
sub-sistema 2, então dE, =(dW)Rev. =(f. ds)=-P,dV, 
em que f é a força exercida na fronteira do sistema e P, 
é a pressão do sub-sistema 1), 


dE=dE, +dE,=-pdV=-P, dV,, (6) 


te 
(dado estarmos a considerar S = c — 


edV=dV,+tdV,=dV), 


temos 
"dE, 9E, 5E, p 
dE = 3, dVi+tgs,. dV, + E dsS,, (7) 
| 9E | E, | 
PF. dVi=75, dS,-P, dVitgs, dS.. (8) 


Como se salientou S= S, + S,e dS = o (reversível). 
Logo 


dS,=-dS, (9) 
e portanto 

E,  0E, 

o, 05, (10) 


Como (E,.S,) são constantes e podemos admitir 
ame de 2) que 


=m E, (11) 

* 4=m58, (12) 
temos 

E,  0E., | 

o “E (13) 


Temos, portanto, que 


E, — te 
5, —€ (14) 


ao longo dos pontos da transformação em que 
o sub-sistema | contacta 2, 


ES = 7 
Designemos esta constante característica de 2 por l 
E,  9E, 


qo “8 


Poderemos, portanto, escrever genêricamente a se- 
guinte relação: 


E=E(S,V), 
—- 9E 9E 
de= dVt5g ás, 
dE=-pdV+LdS. (15) 


A Temperatura Kelvin 
p 


Consideremos um 
ciclo de Carnot 
definido pelos pontos 
PQRS (Fig. 2). 


FIG. 2 


Entre P e Q o sub-sistema que descreve o ciclo con- 
tacta com uma fonte F, caracterizada por £,. 

Entre S e R o sub-sistema contacta F, caracterizada 
por 6. 

As transformações QR e PS são isentrópicas carac- 
terizadas por S, e S, respectivamente. 

Entre P e Q temos (consideramos que É representa 
agora a energia do sub-sistema que descreve o ciclo e 
que entre P e Q contacta F,): 


dE =-pdV + ds, (16) 
dE = - pdV + & dS (contacto com F, 
entre Re sS). (17) 
O rendimento do ciclo define-se pela relação 
— lol (tê) 
| Er | 


em que «w representa o trabalho realizado no ciclo e 
AEr, representa a variação de energia da fonte F,. 


Se chamarmos à energia do sistema E, e S à entro- 
pia, temos que 
E=E+E, +Er, (19) 


em que E é a energia do sub-sistema que descreve o 
ciclo e E,, e Er, são respectivamente as energias das 
Fontes Fe F,(S=S+ S; + Sp). 

No ciclo verifica-se a seguinte relação 


faB=[dE, + JdE, (dadofdE=0) (20) 
PQ RS 


em que admitimos o ciclo a ser descrito no 
sentido PQRSP. 
Logo 
E = ÕE, e AEp, (21) 


(PORSP) PQ RS 


Simplificando a notação temos que 


AE = AE, + 4Er, (22) 
Como dE = - pdV (veja relação (6) ) 
J- pdV = =|- pdV a: dS + SE, (23) 
PQ 
dado que 
| de=] ae + de, (24) 
P Q 


e dE =-pdV+ldS (L—& ao longo de PQ). 
Deste modo 


£, ASpo = - Er, (25) 
e da mesma forma 
[ ASas = - AE. (26) 


Temos, portanto, a seguinte relação: 
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6, ASpo o SEr, | (27) 


[OS AB, 


Como e = [ é maior que zero [2] 
6, |ASpo| ” [AEs | . (28) 
EJASas| [Er 

e como of dS= 0 = ASpo + ASps 

E (29) 

5, [AEr,] 

O módulo do trabalho | «w | é dado por 
dE = dE+dE,. (30) 

dio O .nav=| dE dE... 

b dE b pdV | a /. F», (31) 
w= Er + AE, , (32) 
-[w | =-|4E; | + |4Es |, (33) 


dado que o trabalho é negativo e ÕEr, é negativo e 
ÕEr é positivo no ciclo descrito no sentido horário. 
Como (> 0€e ASpo — - ASas temos que de (27) ÕE,, e 
ÕE,, têm sinais contrários. Se É, > 6, AE, | > AE, E 

Logo verifica-se (33). 


Deste modo 
= | ii | Ra AE, | 
q = [AE | O SE, | ni 
1 
n=l- - dado (29). (35) 
| 


Históricamente a temperatura absoluta Kelvin T defi- 
ne-se pela seguinte relação estabelecida através dum ciclo 
de Carnot em que as fontes têm temperaturas empíricas 
0,e0, 


E 


= 1. = 36 
n T, (36) 
Deste modo podemos fazer 
= 1 (37) 
e passar a escrever 
dE =-pdV+TdsS (38) 


designando T por temperatura Kelvin. 

De notar que a relação (38) é válida numa trans- 
formação infinitesimal (quasiestática) e não necessária- 
mente reversível. 

De facto 


dE = DE gy 


(d - diferencial exterior) 
Ox 


dE =< dE,dP> (contracção da forma dE e 
do vector dP onde dP é dP 
= dx'e ; que pode ter uma 
direcção qualquer e não ser 
necessáriamente tangente à 
transformação reversível). 


Como E = E(S,V) 


dE=(57) AS + (SVO dv. 


De notar também que se E se referir ao sistema 
e a transformação for reversível 


dE = Es dV dadodS=0. 


No entanto para um sub-sistema e numa reversível 
9E 
dE, = A AV, + as ds, 
onde dS, pode ser maior ou menor que zero. 

A variação de entropia é nula numa transformação 
reversível quando nos referimos ao sistema e pode ser 
positiva ou negativa quando nos referimos a um sub-sis- 
tema. 

Dado o anteriormente afirmado podemos enunciar 
com generalidade o Segundo Princípio da seguinte ma- 
neira: 

A entropia dum Sistema só pode aumentar ou per- 
manecer constante (*). 


CONCLUSÃO 


Nos últimos anos tem-se assistido a uma cuidadosa 
revisão conceptual dos Princípios da Termodinâmica. 

Conforme se pode verificar nos artigos referidos 
[3,4,5] o conceito de temperatura surge sem se pôr em 
causa o Primeiro Princípio da Termodinâmica, A revi- 
são conceptual do Primeiro Princípio foi feita pelo autor 
[1]. As conclusões então obtidas são as de que com 
generalidade não é possível atribuir sentido físico à ex- 
pressão de Clausius do Primeiro Princípio dE = dW + 
+ dQ. Sendo assim deve ser possível construir conceptual- 
mente a termodinâmica por outra via diferente da con- 
sagrada por longos anos de tradição. Para tal torna-se 
necessário evidenciar de forma directa a relação entre a 
energia e a entropia [2]. 

Estabelecida esta relação a temperatura surge com 
grande facilidade como um conceito derivado não tendo 
que ser imposta à partida. 

Adquire-se deste modo grande simplicidade e eco- 
nomia de conceitos. 


(*) O enunciado, frequentemente considerado geral, de que a varia- 
ção de entropia de um sistema isolado é maior ou igual a zero é, 
evidentemente, um caso particular [6]. 


- Rodrigo de Abreu, Contribuição para o esclarecimento de al- 


guns problemas fundamentais associados aos primeiro e segundo 
princípios da termodinâmica — transformações relativistas das 
grandezas da termodinâmica, Tese, Instituto Superior Técnico, 
Lisboa (1983) (não publicado). 


?. Rodrigo de Abreu, Variação de Entropia e Reversibilidade, 
| Técnica, N.º 1/2, 1985. 


- D. Home, Concept of temperature without the Zeroth law, 


Am. J. Phys, 45(12), Dec. 1977. 


- P. Ehrlich, The Concept of temperature and its dependence on 


the laws of thermodynamics, Am. J. Phys., 49(7), July 1981. 


« J.S. Thomsen, Comment on “The Concept of temperature and 


its dependence on the laws of thermodynamies”, Am. J. Phys. 
51(5), May 1983. 


. Rodrigo de Abreu, O Conceito de Sistema e o significado físico 


das grandezas trabalho e calor (a ser publicado). 


119 


ON THE VARIANCE OF THE OUTPUT 
SIGNALS FROM SPECTRAL ANALYZERS 


F.E. REBELO SIMÕES (*) 


ABSTRACT 


An exact expression for the normalized variance of 
the output signal from spectral analyzers is derived, for 
the special case of nonheterodyne analyzers with cons- 
tant percentage bandwidth ideal filters and RC integra- 
tor. The input is assumed to consist of large samples of 
a stationary Gaussian process. 


RESUMO 


Determina-se a expressão exacta da variância nor- 
mada dos sinais de saida do analisador espectral não 
heterodínico com filtros ideais de banda de percenta- 
gem constante e integrador RC. Supõe-se que os sinais 
de entrada são amostras suficientemente extensas dum 
processo gaussiano estacionário. 


|. INTRODUCTION 


The aim of this short paper is to determine the exact 
expression of the normalized variance of output signals 
from non-heterodyne spectral analyzers for stationary 
Gaussian input samples (fig. 1). 

The analyzer is supposed to consist of a set of cons- 
tant percentage bandwidth bandpass filters with impulse 


+ FILTER 


* Instituto Superior Técnico and Centro de Análise e Processamento 
de Sinais da Universidade Técnica de Lisboa. 
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response h(t), a squaring device and a RC lowpass fil- 
ter (averaging integrator) of impulse response g(t). 
We have thus 


vít) =[. h(u)x(t -u) dy (1) 
yv() = ví (t) (2) 
with z(t) =). g(v) y (t-v) dv (3) 


g(t)= (I/RC)exp(-t/RC) 1>0 
= () t<O0 
The following simplifying hypotheses are assumed: 
(a) the input sample duration is much larger than 
the time constant RC, 
(b) the power spectrum S, (w) of the input process 
is smooth at filter bandwidth scale. 
(c) each filter is ideal rectangular. 
For the assumed filters the following parameters are 
defined (fig. 2) 


(4) 


- 


| 


FR 
FIG. 2 


F,, F, - lower and upper limit frequencies 
F - center frequency (geometric average) 
F, - medium frequency (arithmetic average) 


Vk =F,/F=EFJ/F, 
B=F, - F, bandwidth 
p= 100 B/F bandwidth percentage 


Two examples are studied: the 1/3 octave and the 
1% filter, for which k =/2, p = 23.156%; and k = 
= 1.01005...., p= 1%, respectively. 

The relative differences (F, - F)/F are respectively 
6.68 x 10” and 1.25 x 10, which means that we may 
take, in the following, F, = F. 


2. DETERMINATION OF THE STATISTICAL 
AVERAGE Ez(t) 


Fig. 3 shows the input signal power spectrum and 
the power transfer function of the filter, centered at f = 
=F=E. 

From egs. (2) to (4), and under hypothesis (a), we get 
Ez(t) = Ey(t) = R/0) (5) 
R,(r) being the autocorrelation function of v(t). 

As R,(r) is the inverse Fourier transform of the power 
spectrum S,(w), given (fig. 4) by 

S (0) = S,(w)| H (w) | (6) 


it follows = 
Ex() =) S,(0)| H(o) [É do/2m (7) 


Fia,3 


Expanding S,(w) in Taylor series in the vicinity of 
f= F,itis possible to obtain [1], [2] for smooth spectra 


Ez(t) = 2B[S,(27F) + D] (8) 
with = 
' + [EC | (9) 
— 24 df? F 


Expression (8) suggests the choice of 
S 27 F) = 2(t)/2B (10) 


as an estimator of S,(27F), so that D becomes the 
bias of the estimator. 


3. DETERMINATION OF THE STATISTICAL 
VARIANCE OF z(t). 
APPROXIMATE EXPRESSIONS 


Let z,(t) be the zero mean signal z(t) - Ez(t). 

The autocorrelation function of this signal is the 
covariance function of z(t), C,(r), of which the Fourier 
transform is the power spectrum S, (w). In the same 
way we define y,(t) and S, (w), the two spectra being 
related by 

Sy (0) = So (0) | G(u) | (11) 
where G(w) is the Fourier transform of g(t). 
From the above we may conclude that 


var rt) =C,(0)=) S,(w) |G(w)|º do/27 (12) 


In order to evaluate S, (w) we express it in terms of 
its inverse Fourier transform, C,(7) which may be writ- 
ten, from (2), as 

Cr) = Ev(vi(t+r) - Ev EvV(t+r) 
= Evi(t)v(t+r) - R$(0) (13) 
in which Ev = 0, because of the filtering action. 

Now, for zero mean Gaussian variables one has the 
known relation [3] 

Ev, V;V; Vs= Ev, v; Ev; vs + Ev, v; Ev, vs + Ev, vs Evo v; (14) 
which vields the result 

C,(7) = 2RV(7) (15) 
In the frequency domain and since S(w) is an even 
function, the above relation takes the form 


S(u)=2) SANSAN+Ha)dA/2m (16) 


For a rectangular filter of center frequency F and 
symmetrical about F, in linear scale, the spectrum Sy A) 
has the form represented in fig. 5, under the hypothesis 
(b). It must be remarked that a variation of 35% of 
S,(w) within the filterband will produce only a 1% 
error on, (w). Fig. 5 shows also the graphic of | G(w) |. 

As for the integrator a “noise bandwidth” is defined 


B,=). |G(w)|*do/27 | [G(w) [ma = 1/T (17) 
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LS yolw) 


4BS; (277F) 


FIG.5 


in which T'is the “averaging time”. 

As it will be seen later, it is desirable to have B> B,, 
1.e., BT>H1, for reducing the measuring error. In these 
circumstances, it is usual practice [2] in the evaluation of 
var z(t) to take S, (0) for S, (w) in the integral (12), 
because of the concentiation of |G(w)|” near zero, 
thus disregarding the slope of the central triangle and 
the presence of the other two. 

In this way we get 
var z(t)=4BB,S o7F)=4BS,(27F)/T (18) 


Under hypothesis (b), D < S(27F), so that we can 
write in normalized form 


2 — var z(t) = A 
[Ex)f BI 


which also holds for the estimator S.(2 mF), by (10). 


(19) 


Tº 


1 T OCTAVE 


DRUnaA Pr HE 


o, 
sal 
1Q 


+ RE eh Zi 


BT 
1000 
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2BS x (271) 


4. EXACT EXPRESSIONS OF varz(t) ANDo/[1] 


the cases for which BT is not much greater than unity, 
we only have to substitute in (12) the expression of 
S,w), easily taken from fig. 5, and |G(w)|?, given by 


[G(w)|* = 1/[1+ (vRC)] 
which is readily obtained from (4). 
This last expression also yields the averaging time 


In order to correct the expressions (18) and (19) in 
(20) 


T=2RC (21) 
The final result of the integration is the following: 


var z(t)=4B'S/(27F)[U(BT7m)-V(BTm)] (22) 
o” = U(BT7)- V(BT7) (23) 


BTO* 


TIUOIAO 
Damian 


RP O O 


| EE 
ES 


ni 
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NICE 
PISSACOTO à 


100 1000 | 


00001 0.001 001 
FIG 7 


where 

wU(w) = 2arctg w - 2m arctg mw + (m+1) arctg 
(m+lw + (m-l)aretg (m-l)w 

wV(w) = log(1+w?) - log[lH-m?w?] + 

Alog [H(m+5 log [H+(m-1)'w?] (24) 


m=2F,/B (25) 


We can conclude from (23)/(24) that 


I BTo'= 

BT (26) 
lim o” = 

BT>O À am 


Figs. 6 and 7 show o” and BTo” as functions of BT, 
for the two filters of 1/3 octave and 1%, for which m = 
8.6946... and 200,0025... respectively. For an easier exam- 
ination of the curves a logarithmic scale was used. 

We may conclude that in either case it is only legi- 
timate to use the approximation 1/BT for relatively 
large values of BT. 

In fact, for BT = 10 the approximation 1/BT yields 


an error of 9%, while for BT = 5 the error rises to 15%. 
If we go down further we find that the true value of o” 
is 0,56 instead of |, for BT = 1. 

For extremely low values of BT, the above appro- 
ximation becomes meaningless. In fact, for BT = 0.05 
we will get for o” the true values 1.346 and 0.997 for 
1/3 octave and 1% filters, against the approximation 
value 1/BT = 20. 

As a final comment we must note that in both appro- 
ximate and exact expressions, the admissible values of 
o” for a reasonable measurement imply relatively large 
values of BT, as we stated before. 
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NORMAS PARA ADMISSÃO 


DE TRABALHOS E 


INSTRUÇÕES PARA OS AUTORES 


1. A Revista TÉCNICA, em princípio, publicará anualmente 
quatro números de natureza técnico-científica. 


2. Atendendo ao espaço disponível e carácter científico que se 
deseja imprimir à Revista, as três formas de publicação previstas, 
revestir-se-ão das seguintes características: 


ARTIGOS POR CONVITE — Trabalho de sintese, no domínio 
da investigação científica ou no domínio do desenvolvimento, 
contemplando actividade desenvolvida entre nós, ou pelo menos, 
de interesse para o desenvolvimento científico e tecnológico do 
Pais. 

ARTIGOS TÉCNICO-CIENTÍFICOS — Trabalho de investiga- 
ção científica ou de aplicação tecnológica a um domínio bem 
definido da ciência ou da técnica e de que o autor seja um especia- 
lista. 

NOTAS CIENTÍFICAS — Destinam-se essencialmente à apre- 
sentação de novos resultados no domínio da investigação cienti- 
fica ou do desenvolvimento tecnológico. Pretende-se uma nota 
curta, realçando os resultados alcançados e referenciando conve- 
nientemente trabalhos afins, evitando a exposição longa de ante- 
cedentes, que deverá ser remetida para referências. 


3, Os artigos por convite e os artigos técnico-cientificos, deverão 
ser escritos em língua portuguesa. As notas científicas poderão ser 
escritas em língua estrangeira, de preferência em inglês. 


4. Os autores deverão atender na elaboração dos artigos às se- 
guintes condições: 

— os artigos convidados não poderão exceder 30 páginas dacti- 
lografadas, incluindo o espaço reservado a ilustrações; 

— os artigos técnico-científicos não poderão exceder 15 páginas 
dactilografadas, incluindo o espaço reservado a ilustrações; 

— as notas científicas não poderão exceder 8 páginas dactilogra- 
fadas. 


5. Os trabalhos devem ser dactilografados a dois espaços, em 
folhas modelo A4, devidamente numeradas. 


6. Os trabalhos devem ser enviados em duplicado, recomendan- 
do-se aos autores que conservem uma cópia. 
7. Conjuntamente com o trabalho, devem ser enviados: 


— título do trabalho em português e inglês; 

— resumo e abstract; 

— nome, direcção e telefone do autor e demais dados que facili- 
tem o processo de revisão de provas. 
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8. Os autores deverão assegurar-se, antes de enviarem os origi- 
nais, que estes não contenham erros, falta de dados, etc., pois que, 
na fase de revisão de provas, apenas se corrigirão erros de compo- 
sição tipográfica. 


9. As referências bibliográficas serão indicadas no texto por or- 
dem numérica, usando algarismos árabes entre colchetes, e a lista 
das mesmas, no final do artigo, obedecerá às seguintes normas: 
Artigos de revistas: nome(es) do autor(es) e iniciais nome da re- 
vista, número do volume, página e ano, entre parêntesis. O título 
do artigo só deverá ser mencionado se isso constituir matéria 
fundamental para esclarecimento. 

Livros: nome do autor e iniciais, título, editora, ano de publica- 
ção e páginas. 


10. Os Quadros ou Tabelas, numerados em romano, devem figu- 
rar no texto datilografado, podendo ser encimados por um título 
conciso. 


H. Os gráficos devem ser desenhados a tinta da china, em folhas 
separadas e prontos para serem reproduzidos. Quaisquer letras 
ou algarismos que neles figurem, deverão ser suficientemente 
grandes para permitir a sua redução que, em princípio, fica ao 
critério do editor. 

As folhas com os gráficos, figuras ou fotografias, serão numera- 
das, indicando-se o trabalho a que dizem respeito. Além disso, no 
texto, indicar-se-à o local para a respectiva inclusão. 


NOTA: 4 TÉCNICA salienta a necessidade de os autores, 
quando for o caso, obterem em tempo as devidas autorizações de 
reprodução de excertos de outras obras, sendo esta matéria da 
inteira responsabilidade do autor. 


A TÉCNICA procurará, sempre que possível e por princípio, o 
parecer de individualidades de reconhecida competência técnico- 
científica, no sentido de apurar da originalidade, qualidade e inte- 
resse científico dos trabalhos que lhe são enviados. Esta atitude 
não se aplica aos artigos convidados, a menos que os autores o 
desejem. À Direcção da TÉCNICA, caberá sempre a ceitação ou 
recusa dos originais que lhe sejam enviados. A TÉCNICA notifi- 
cará os autores ou devolverá os trabalhos submetidos a aprecia- 
ção num prazo máximo de 90 dias após a recepção do original. Os 
trabalhos publicados, serão registados com a data de recepção do 
original na TÉCNICA e disso se dará pública notícia em respec- 
tiva nota de fim de página. 

A TÉCNICA oferecerá cinco exemplares da Revista, por cada 
trabalho publicado, independentemente do número de autores. 


DO a q | 


anuncie 
na técnica 


Com uma tiragem de 5 000 exemplares e uma difusão 

à escala nacional e internacional a “TÉCNICA” abarca nos seus leitores 
quadros científicos e técnicos das mais variadas empresas, bem comc 
grande número de professores e estudantes universitários. 

Face a isto, poderemos concluir que se a “TÉCNICA” 

merece o apoio de todos os que de algum modo se encontram 
ligados à Engenharia — Escolas, Institutos de Investigação e Empresas 
— ela é, para além de um veículo da investigação teórica e experimental, 
um veículo por excelência de produtos e serviços, em suma, 

das actividas da vossa empresa. 

Através da inserção de publicidade da vossa empresa na revista 
“TÉCNICA” está apoiando esta publicação e comunicando 

com um universo de potenciais interessados nas vossas actividades. 


ANUNCIE NA “TÉCNICA”, CONTACTE-NOS: 
Técnica . aeist 
Av. Rovisco Pais — 1000 LISBOA — Tels. 88 10 18/8975 15 
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